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Abstract. We study here a number of mathematical problems related to 
our recently introduced neoclassical theory for electromagnetic phenomena in 
which charges are represented by complex valued wave functions as in the 
Schrodinger wave mechanics. In the non-rclativistic case the dynamics of el- 
ementary charges is governed by a system of nonlinear Schrodinger equations 
coupled with the electromagnetic fields, and we prove that if the wave functions 
of charges are well separated and localized their centers converge to trajecto- 
ries of the classical point charges governed by the Newton's equations with the 
Lorcntz forces. We also found exact solutions in the form of localized acceler- 
ating solitons. Our studies of a class of time multiharmonic solutions of the 
same field equations show that they satisfy Planck-Einstein relation and that 
the energy levels of the nonlinear eigenvalue problem for the hydrogen atom 
converge to the well-known energy levels of the linear Schrodinger operator 
when the free charge size is much larger than the Bohr radius. 



1. Introduction 

It is well known that the concept of a point charge interacting with the electro- 
magnetic (EM) field has fundamental problems. Indeed, in the classical electrody- 
namics the evolution of a point charge q of a mass m in an external electromagnetic 
(EM) field is governed by Newton's equation 



^ [mv (t)] = q 



E(t,r(t)) + -v(f) xB(i,r(t)) 
c 



(1) 



where r and v = r =4| are respectively the charge's position and velocity, E(i,r) 
and B (t, r) are the electric field and the magnetic induction, and the right-hand 
side of the equation (1) is the Lorentz force. On other hand, if the charge's time- 
dependent position and velocity are respectively r and v then the corresponding 
EM field is described by the Maxwell equations 

1 $B 

-— + V x E = 0, V ■ B = 0, (2) 

c ot 

~-VxB = --?*(x-r(t))v(t), V-E = 47r«*(x-r(t)) ) (3) 
c ot c 

where S is the Dirac delta- function, v (t) — r (t), c is speed of light. If one would 

like to consider equations (l)-(3) as a closed system "chargc-EM field" there is 

a problem. Its origin is a singularity of the EM field exactly at the position of 
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the point charge, as, for instance, for the electrostatic field E with the Coulomb's 



electromagentic theory, a possible remedy is the introduction of an extended charge 
which, though very small, is not a point. There are two most well known models 
for such an extended charge: the semi-relativistic Abraham rigid charge model (a 
rigid sphere with spherically symmetric charge distribution), [35, Sections 2.4, 4.1, 
10.2, 13], [32, Sections 2.2], and the Lorentz relativistically covariant model which 
was studied and advanced in [1], [20, Sections 16], [28], [30], [32, Sections 2, 6], [33], 
[35, Sections 2.5, 4.2, 10.1], [37]. Poincare suggested in 1905-1906, [31] (see also [20, 
Sections 16.4-16.6], [32, Sections 2.3, 6.1- 6.3], [29, Section 63], [33], [37, Section 4.2] 
and references there in) , to add to the Lorentz- Abraham model non-electromagnetic 
cohesive forces which balance the charge internal repulsive electromagnetic forces 
and remarkably restore also the covariance of the entire model. 

Another well known problem of a physical nature with point charges and the clas- 
sical electrodynamics is related to the Rutherford planetary model for the hydrogen 
atom. The planetary model of the hydrogen atom is inconsistent with physically 
observed phenomena such as the stability of atoms and the discreteness of the hy- 
drogen energy levels. That inconsistency lead as we well know to introduction of 
non-classical models, such as the Bohr model and later on the Schrodingcr's model 
of the hydrogen atom. 

To address the above mentioned problems with the classical electrodynamics we 
introduced recently in [5], [6] wave-corpuscle mechanics (WCM) for charges which 
is conceived as one mechanics for macroscopic and atomic scales. This model is 
neoclassical in the sense that it is based on the classical concept of the electromag- 
netic field but an elementary charge is not a point but a wave-corpuscle described 
by a complex- valued scalar wave function ip with the density \ip (i,x)| 2 which is not 
given a probabilistic interpretation. The dynamics and shape of a wave corpuscle 
is governed by a nonlinear Klein-Gordon or a nonlinear Schrodinger equation in 
relativistic and nonrelativistic cases respectively, and a wave-corpuscle is defined as 
a special type of solutions to these equations. Consequently an elementary charge 
in the WCM does not have a prescribed geometry which differs it from the classical 
Abraham and Lorentz models and more recent developments (see [7], [23], [25], [18], 
[35] and references therein). Another important difference is that Newton's equa- 
tions are not postulated as in (1) but rather arc derived from the field equations in 
a non-relativistic regime when charges are well separated and localized. Physical 
aspects of the model are discussed in detail in [5], [6], where the Poincare type 
cohesive forces are constructed somewhat differently then we propose here. The 
primary focus of this paper is on studies of mathematical properties of the field 
equations of the WCM. 

The paper is organized as follows. In the first two subsections of this Intro- 
duction we describe the system relativistic and non-relativistic Lagrangians and 
the corresponding field equations, and in the third one we consider the relevant 
nonlincarities and their examples. Importantly, the nonlinearity depends on a size 
parameter a > associated with a free solution (ground state) . Section 2 is devoted 
to remote interaction regimes for many charges. We prove there that in the case of 
the non-relativistic field equations when a — > (macroscopic limit) the centers of 
the interacting charges defined by the formula 



potential 



with a singularity at x = r. If one wants to stay within the classical 
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converge to the solutions of Newton's equations with the Lorentz forces if ip remain 
localized. We also provide examples of exact solutions of the field equations in 
the form of accelerating solitons for which the localization assumption holds. In 
Section 3 we study systems of bound charges and time multiharmonic solutions to 
the field equations. We consider there, in particular, a connection, discovered in 
a different setting in [12], between the Planck-Einstein energy- frequency relation 
and the logarithmic nonlinearity as well as some dynamical issues related to the 
logarithmic nonlinearity closely related to results in [14], [16]. We continue then 
with a study of a system of two charges with the logarithmic nonlinearity as a model 
for the hydrogen atom. We prove that if n = a±/a — > 0, where a\ is the Bohr radius, 
then the lower energy levels associated with this model converge to the well known 
energy levels of the linear Schrodinger operator for the hydrogen atom. The proof 
is based on an approach developed in [9], [10] with certain modifications. 



1.1. Relativistic Lagrangian and field equations. Let us consider a system 
of N charges interacting directly only with the EM field described by its 4-vector 
potential (<p,A.). The charges are described by their wave functions ip e with the 
superscript index £ = 1,...,N labeling them. The EM fields are related to the 
potentials 92, A by the following standard relations 



-V<p - -St A, B = V x A. 

c 



(4) 



We introduce now the system Lagrangian C which is though similar to the one in- 
troduced in [5], [6] but differs from it. A different Lagrangian is introduced here in 
order to get the hydrogen atom frequency spectrum with desired precision. Namely, 
to every £-th charge is assigned an adjunct potential (<//,A^) describing its addi- 
tional degrees of freedom and our relativistic Lagrangian is defined by 



C 



Vip+^dtA) -(VxA) 2 



(5) 



E 



X 
2m l 



2 r 



-2 r* 



E- 



K 2 



G' 



(W) 



where (i) the covariant derivatives are defined by 



iq e (A - A') 



X c 



(6) 



(ii) tp* is complex conjugate to ip; (iii) m > is the £-th charge mass, q e is the value 
of the charge, kqi = S^ 2 , and x > is a constant similar to the Planck constant 

H = (iv) G l (ip l *ip e ) is a nonlinearity which will be described below. 
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The Euler-Lagrange field equations for the Lagrangian C defined in (5) include, 
first of all, the Maxwell equations for the EM potentials 



d t A + V<£ 



V x (V x A) + -d t ( -dtA + V<p) = -Vj 



47T 



where the charge densities and currents are defined by 

„ 2 



nv 



Xlm lF c 



(7) 

(8) 

(9) 
(10) 



The field equations include similar Maxwell equations for the adjunct potentials 

1 



-d t A + Vip' 



1 



V x (V x A 1 ) + -d t -d t A e + Vtp 



1 



-4irp\ 

477 



-y, e = i,...,N, 



(ii) 

(12) 



and equations for the wave functions ip in the form of nonlinear Klein-Gordon 
equations 

2~' 



fdfdW + V V - G" (VV ) V £ - ^ l = 0, l = 1, N. (13) 



Note that equations (13) for ip are coupled to the equations for EM potentials via 
the covariant derivatives. 

The Lagrangian and equations are Lorentz and gauge invariant, and one can 
verify that defined above p and J satisfy conservation/continuity equation 

-8^ + v • j' = o. 

c 

If we choose the Lorentz gauge 

-dtp + V • A = 0, -d t w £ + V • A e = 
c c 

for all potentials, equations (7)-(12) take the form 

n 



v^-4stV = -EW, 
c 1=1 

ld 2 A-V 2 A = ^fy, 



(14) 

(15) 
(16) 



vv - 4 a tV = -w, 

c z 

la t 2 A^-V 2 A £ = — J i ,£ = l,...,N, 
c z c 



(17) 
(18) 
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where p l , 3 l are defined by (9), (10). Based on the above equations we assume that 
always 



(<p,A) = Yl<(fi,A<). (19) 



Hence, from now on we assume that in the relativistic case the dynamics of EM fields 
and charges is determined by equations (17), (18), (9), (10), (13) and according to 

(19) and (6) the covariant derivatives d\ and in (13) involve only fields (w e ' , A e '^j 

with (! ^ I, 



# = * + v' = V-^- (20) 

Y YC 

implying that the adjunct potentials completely compensate the EM self-action for 
every charge and effectively there is no EM self- interaction (in [5], [6] this kind of 
compensation was gained by using an additional nonlinear self-interaction). 

Let us turn now to the nonlinearities G. We introduce the nonlinearities here 
the same way as in [5], [6] with an intension to have a rest (ground) state for a 
single charge in the form of a localized wave function (form factor). Notice that for 
a single charge evidently N = 1, <p — <// = 0, A — A e = and we look for the rest 
solution in the form 

(t, x) = e- iw °V (x) , loo = — = ck (21) 

X 

with A e = 0. Substituting (21) in (13) yields 

-VV = 4tt|V>| 2 , 

-V 2 ip + G' (V>»V = °- 

Now we choose a strictly positive, monotonically decreasing radial function ip (ground 
state) as a parameter for the model and determine the nonlinearity G 1 from the fol- 
lowing charge equilibrium condition: 



V 2 V> + G' ' 



i\) = 0. (22) 



The above equation allows to determine the nonlinearity G' as long as tp is a strictly 
positive, smooth and monotonically decreasing function of r = |x|. In Section 1.3 
we consider the nonlinearity in more details providing also examples. 

Note that using the Lorentz invariance of the system one can easily obtain a 
solution which represents the charge-field moving with a constant velocity v simply 
by applying to the rest solution ( ip e , (p,0j the Lorentz transformation (see [5], [6]). 



1.2. Non-relativistic Lagrangian and field equations. Our non-relativistic 
model describes the case of charges moving with non-relativistic velocities and it 
is set as follows. Using the frequency-shifting substitution (21) with more general 
if) = ip^ (i, x) which depends on (t, x) we observe that the second time derivative in 
(13) can be written in the form 



(i 
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where 



(23) 



We neglect the term with the factor 4j and substitute ~2i^- [d t + ^-^P^e^j ^Pu + 

KQipu for the term —^d^dfiplj in (13). Consequently, we replace the nonlinear Klein- 
Gordon equation (13) by the following nonlinear Schrodingcr equation (where we 
denote the frequency shifted ip^ by ip ) 



X 

2m e 



2m e 



G 



where V is given by (20). Since the magnetic fields generated by moving charges 
also have coefficient - we neglect them preserving only the external magnetic fields 

- £ ~ £ 

and replace V by the covariant gradient V ox defined by 



So our non-relativistic Lagrangian is 

N 



ft - 



|vg|- 

8tt 



(25) 



(26) 



£ I i £\ i£ /£* 

-q W + Vex - f ) ^ i> 



8tt 



where A ex (t, x) and ip cx (t, x) are potentials of external EM fields, ip * is complex 
conjugate of ip £ . The Euler-Lagrange equations for the electrostatic potentials have 
the form 

ivV = f>W*, (27) 



1 



-vy = 9 Wv = i, 



iv. 



(28) 



Assuming that that <^ vanish at infinity we obtain a reduced version of (19) 

¥> = Xy. (29) 

Similarly to the relativistic case we assume now that nonrelativistic equations for 
dynamics of charges in the external EM field with potentials </? ex , A ex take the form 

= -£- t (vL) V + 9 £ + v cx ) ^ + ^ [<*]' (|^f) ^, (30) 



where is given by (23) and <p is determined from the equations 

2 



v V = -W 



= 1,...,7V. 



(31) 
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The solution of the above equation is given by the formula 

2 



<f (t,x) = <f I -dy. (32) 



(t,y) 



x 



The nonlinear self-interaction terms G a in (30) are determined through the charge 
equilibrium equation (22) and index a indicates the dependence on the size param- 
eter a > which we introduce by the formula 

G' a (s) = a- 2 G[ (a 3 s) . (33) 

More detailed discussion of the nonlinearity is given in the following section. 

1.3. Nonlinearity, its basic properties and examples. As wc have already 
mentioned, the nonlinear self interaction function G is determined from the charge 
equilibrium equation (22) based on the form factor (free ground state) ip. Important 
features of our nonlinearity include: (i) the boundedness or slow subcritical growth 
of its derivative G' (s) for s — > oo with consequent boundedness from below of the 
energy; (ii) slightly singular behavior about s = 0, that is for small wave amplitudes. 

In this section we consider the construction of the function G, study its properties 
and provide examples for which the construction of G is carried out explicitly. 
Throughout this section we have 

ip,ip > and hence \%p\ = tp. 

We introduce explicitly the dependence of the free ground state ip on the size pa- 
rameter a > through the following representation of the function ip (r) 

^W = iW^" 3/2 M 8 " 1 ''). ( 34 ) 
where i/> x (r) is a function of the dimcnsionless variable r > 0. The dependence on 



a is chosen so that L 2 -norm 



docs not depend on a, hence the function 



■0 a (r) satisfies the charge normalization condition (53) for every a > 0. Obviously, 
definition (34) is consistent with (22) and (33). The size parameter a naturally has 
the dimension of length. A properly defined spatial size of tp a , based, for instance, 
on the variance, is proportional to a with a coefficient depending on ip ± . The charge 
equilibrium equation (22) can be written in the following form: 

V 2 ^ a = G' a ], a . (35) 

The function ip a (r) is assumed to be a smooth (at least twice continuously dif- 
fcrentiablc) positive monotonically decreasing function of r > which is square 

integrable with weight r 2 , we assume that its derivative "0 a (r) is negative for r > 
and we assume it to satisfy the charge normalization condition of the form (53); 
such a function is usually called in literature a ground state. 

Let us look first at the case a = 1, i\) a = ipi, tp a = $>i, for which the equation 

o 2 

(35) yields the following representation for G'^i) from (35) 

G ;(^ (r) ) = (Z^Hr). (36) 
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Since •ip 1 (r) is a monotonia function, we can find its inverse r = r (ip 2 ) ; yielding 

G\ (s) = = # (oo) < s < # (0) . (37) 

(»"(*)) 

Since (r) is smooth and d r il,\ < 0, G'(|^>| 2 ) is smooth for < \ip\ 2 < (0). If 

o 2 

we do not need G' (s) to be smooth, we extend G' (s) for s > tp 1 (0) as a constant, 
namely 

G\ (s) = G\ (0)) i£s>i>l (0) • (38) 

° 2 

The first derivative of such an extension at s = ^ (0) has a discontinuity point. If 
tp a (r) is a smooth function of class G™, n > 2, we always can define an extension 

Q 2 

of G' (s) for s > ip 1 (0) as a bounded function of class G™~ 2 for all r > and 

G; (s) = G; (0)) - 1 if a > i,\ (0) + 1. (39) 

Slowly growing (subcritical) functions G' (s) which are not constant for large s also 
can be used, see examples below. 

In the case of arbitrary size parameter a > we define G' a (s) by formula (33), 
and this definition is consistent with (34) and (37). 

Let us take a look at general properties of G' (s) as they follow from defining 
them relations (37). In the examples below the function G' (s) is not differentiable 
at s = 0, but if tp (r) decays exponentially or with a power law the nonlinearity 
g (ip) = G'(|?/>| )ip as it enters the field equation (30) is differentiable for all ip 
including zero, hence it satisfies the Lipschitz condition. For a Gaussian ip 1 (r) 
which decays superexponentially G'(|V>| 2 ) is unbounded at zero and g (ip) is not 
differentiable at zero. Since ^ (|x|) > 0, the sign of G[ (j^l 2 ) coincides with the 
sign of V 2 ^ (|x|). At the origin x = the function ip 1 (|x|) has its maximum and, 

° 2 

consequently, G[ (s) < for s close to s = ip 1 (0). The Laplacian applied to the 
radial function ^ takes the form i £ (r& |x|) . Consequently, if r& (r) is convex 

at r = |x| we have V 2 i/j 1 (|x|) > 0. Since r 2 %p 1 (r) is integrable, we naturally assume 
that |x| ip 1 (|x|) — > as |x| — > oo. Then if the second derivative of rip 1 (r) has a 
constant sign near infinity, it must be non- negative as well as G^ (s) for s< 1. In 
the examples we give below G^ (s) has exactly one zero on the half-axis. 

Example 1. Consider a form factor (r) decaying as a power law, namely 

<Mr) = -2— — =-tj , (40) 

11 ' (l + r 2 ) 5/4 1 ' 

where c pw is the normalization factor, c pw = 3 1 / 2 / (47r) 1 ^ 2 . This function evidently 
is positive and monotonically decreasing. Let us find now G' (s) based on the 
relations (37). An elementary computation of V 2 ^ shows that 

15s 2 / 5 45s 4 / 5 
G(S) = 177T-^' ( 41 ) 

^OLpw T:t-pW 
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The extension for s > c 2 w can be defined as a constant or the same formula (41) 
can be used for all s > since corresponding nonlinearity has a subcritical growth. 
If we explicitly introduce size parameter a into the form factor using (34), we 

o 2 

define G' a (s) by (33). Notice that the variance of the form factor ip 1 (|x|) decaying 
as a power law (40) is infinite. 

Example 2. Now wc consider an exponentially decaying form factor tp 1 of the 
form 

^ 1/2 

fa (r) = c c e-( r2+1 ) 1/2 , c c = UttJ^ rV 2 ^ 1 ) 1 " drj . (42) 

Evidently tp 1 (r) is positive and monotonically decreasing. The dependence r (s) 
defined by the relation (42) is as follows: 

r = [In 2 (cA/i) - 1] V2 . if v 7 ^ < i>! (0) = c c e'\ (43) 
An elementary computation shows that 

v2 ^i 2 1 1 

T 1 = r + m~T\ + ~ - L 44 

iP 1 (r 2 + l)5 (r 2 + l) ( r 2 + 1 )3 

Combining (43) with (44) we readily obtain the following function for s < c 2 e -2 

4 4 8 



G[ (s) 



1 



(45) 



In(cg/a) ln 2 ( C 2/ s ) In 3 (cg/«). 

Wc can extend it for larger s as follows: 

G' 1 ( S ) = G' 1 (^e- 2 )=-3if S >c^e- 2 , (46) 

or we can use a smooth extension as in (39). The function G[ (s) is not differentiable 
at s = 0. At the same time the function g (ip) = G[ {ip (r)) ip if we set g (0) = 
is continuous and <? (ip) is continuously differentiable with respect to ip at zero and 
g ("0) satisfies a Lipschitz condition. The variance of the exponential form factor 
ip l (r) is obviously finite. To find G' a (s) for arbitrary a we use its representation 
(33). 

Example 3. Now wc define a Gaussian form factor by the formula 

^(r) = C g e- r2 /*, C g = -^. (47) 
Such a ground state is called gausson in [12]. Elementary computation shows that 
V 2 V> (r) 



r 2 -3 = -ln(^ 2 (r)/C 2 ) -3. 
ip(r) v / 

Hence, we define the nonlinearity by the formula 

G'(|V1 2 )=-ln(|V1 2 /C 2 )-3, (48) 

and refer to it as logarithmic nonlinearity. The nonlinear potential function has the 
form 

G(s) = jf (- In (s'/C 2 ) - 3) ds' = -s In s + s ^ln - 2^ . (49) 
Dependence on the size parameter a > is given by the formula 

G' a (|V| 2 ) = In (a 3 |^| 2 /C 2 ) - Sa" 2 . (50) 



10 



ANATOLI BABIN AND ALEXANDER FIGOTIN 



Obviously g (iji) = G[ (\ip\ 2 ) ij) is continuous for all ifi £ C if at zero we set g (0) = 
and is differentiable for every ip ^ but is not differentiable at ip = and does not 
satisfy the Lipschitz condition. Notice also that g (ip) has a subcritical growth as 
— > oo. 



2. Charges in remote interaction regimes 

The primary focus of this section is to show that if the size parameter a — > 
then the dynamics of centers of localized solutions is approximated by the Newton 
equations with the Lorcntz forces. This is done in the spirit of the well known in 
the quantum mechanics Ehrenfest Theorem, [34, Sections 7, 23]. We also provide 
as an example explicit wave-corpuscle solutions which have such a dynamics. 

As a first step we describe basic properties of the classical solutions of (30), (31). 
The Lagrangian £ is gauge invariant and every ^-th charge has a 4-current (p e ,J e ) 
defined by 



p = q 



0< 



— t im ■ 



■0 £ m t c 

which satisfies the continuity equations dtp 1 + V • J e = < 



4 



i> 1 



mfc 



= 0. 



(51) 



(52) 



Note that J e defined by (51) agrees with the definition of the current (10) in the 
Maxwell equations. Equations (52) can be obtained by multiplying (30) by ip * 
and taking the imaginary part. Integrating the continuity equation we see that 



const and we impose the following normalization condition: 



dx = 1, t > 0, 1 = 1, 



,N. 



(53) 



The motivation for this particular normalization is that this normalization allows 



to converge to a delta-function and, in addition, for this normalization 



the solution to (31) given by the formula (32) converges to the Coulomb's potential 



with the value q e of the charge if 



converges to a delta function. The momentum 



density P for the Lagrangian Cq in (26) is defined by the formula 



(VvV 



Note that so defined momentum density P e is related with the current J in (51) 
by the formula 



P* (t,x) = - F J^ (t,x). 
We introduce the total individual momenta P e for £-th charge by 



P'dx, 



and obtain the following equations for the total individual momenta 



dP^ 

dt 



E c 



1 



x B r 



dx. 



(54) 



(55) 



(56) 
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1 



(57) 



where 

V (t,x) = i I P^ (t,x) = ^J« (t,x). 

The external EM fields E ox , B cx in (56) corresponding to the potentials ip ex , A ox 
are determined by standard formulas (4). Derivation of (56) is rather elementary. 
Indeed, in the simplest case where A ex = we multiply (30) by V?//*, take the real 
part and integrate the result over the entire space using integration by parts. To 

obtain (56) in more involved general case one can similarly multiply (30) by V ip * 
and then integrate the result by parts using some vector algebra manipulation. 

2.1. Newtonian mechanics as an approximation. Let us show now that if the 
size parameter a is small compared to the macroscopic scale of variation of EM fields, 
the charge evolution can be described approximately by the Newton equations with 
the Lorentz forces similar to (1). In this subsection we give a heuristic derivation. 
Conditions under which this kind of derivation is justified arc given in the next 
section. 

We introduce the £-th charge position r e (t) and velocity v l (t) as the following 
spatial averages 



(t) = r e a (t) 



dx, 



(t) 



3 e (t, x) dx, 



(58) 



where the current density 3 l is defined by (51). We show below that a combination 
of the continuity equation (52) with the momentum evolution equations (56) imply 
the following remarkable property: the positions v l (t) satisfy with a high accuracy 
Newton's equations of motion for the system of N point charges if the size parameter 
a is small enough. 

Multiplying continuity equation (52) by x and integrating we find the following 
identities 



dr^ (t) 
dt 



x<9 t 



dx = — 



1 



xV • J 1 dx = 



1 



K 3 



J f dx = v* (t) , (59) 



showing that the positions and velocities defined by formulas (58) are related exactly 
as in the point charge mechanics. Then integrating (54) we obtain the following 
kinematic representation for the total momentum 



P*( t ) = ™ f j«( f x ) dx = mV(t). 
Q Jr3 



(60) 



which also is exactly the same as for the point charges mechanics. Relations (59) 
and (60) yield 



(t) 



(t) 



dP £ 
~~dT 



(61) 



dH dt 

and we obtain from (56) the following system of equations of motion for TV charges: 



dV (t) 
dH 



1 



E f '+E ox ) ib* +-v £ xB cx dx, £ = 1,...,N, 
/ c 

(62) 

where E f (t, x) = — V<p f (t,x), E cx and B cx are defined by (4). 

The derivation of the above system is analogous to the well known in quantum 
mechanics Ehrenfest Theorem, [34, Sections 7, 23], [12]. Now we give a formal 
derivation of Newton's law of motion. Let us suppose that for every ^-th charge 

2 



density 



J 



and the corresponding current density J are localized in a- vicinity of 
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the position v l (t) , and that r e (t) — v l (t) 
time interval [0,T]. Then if a — >• we get 

2 



> 7 > with 7 independent on a on 



(t, x)^8( X -r e (t)) , (t, x) = 3 e /q e / (t) S(x-r e (t)) , (63) 



where the coefficients before the Dirac delta-functions are determined by the charge 
normalization conditions (53) and relations (58). Using potential representations 
(32) we infer from (63) the convergence of the potentials <p to the corresponding 
Coulomb's potentials, namely 

/ (t, x) -> p§ (f , x) = -^--^ , V V (t, x) -> g ^ X ~ y as a -> 0. (64) 
l x — r I |x — r*| 

Hence, when passing to the limit as a — > we can recast the equations of motion 
(62) as the system 



dt 2 



co, (65) 



where 



f = £ «' + g*E ex (/) + ^v^ x B cx (/) , ^ = 1, ...,7V, 



and eo — > as a — > 0. Notice that the terms f f in equations (65) coincide with the 
Lorentz forces and we see that the limit equations of motion obtained from (65) 
coincide with Newton's equations of motion for point charges interacting via the 
Coulomb forces and with the external EM field via corresponding Lorentz forces, 
namely 

l5" = -E | „3 ) + ^ + -/ >< Bex (r 1 ) ,t = l,...,N. (66) 



The above formal derivation of Newton's equations of motion is very simple, but 
if implemented with full rigor would require to impose rather high regularity con- 
ditions on the solution and rather strong requirements on the convergence to the 
delta function. Note that we use essentially the fact that the nonlinearity G' a , that 
singularly depends on a as a — > according to (33), does not enter the system (62). 
In the following subsection we rewrite the system (62) in an integral form which 
allows to pass to the limit under less restrictive assumptions, and also study an 
important case of a single charge in an external field. 

2.2. Single charge in an external EM field. In this section we describe assump- 
tions under which the above derivation of Newton's law with the Lorentz forces for 
localized charges can be rigorously justified. In the case of many well separated 
charges every single charge senses all other charges through the sum of their EM 
field, and it is this feature of the EM interaction makes a problem of a single charge 
in external EM field of special importantce. 

We consider a single charge described by (30) with N = 1 in an external EM 
field where tp cx (t,x), A ex (t,x) are two times continuously differcntiablc functions. 
Now index £ takes only one value, but we still keep it to be consistent with our 
previous notation. We assume that the first spatial derivatives of E ox and B GX are 
bounded uniformly in (t,x). We also assume that equation (62) holds in a weaker 
sense which we describe below. 
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First let us recast (62) in a different form. Multiplying (52) by x — r e (t) we 
obtain 



m l d t 



(x - / (*)) + 



m 



i> 1 



+V 







It 



d t r* (t) + 
(x-r f (t)) 



ylm 



</■' 



Note that the right-hand side of the above equation coincides with m e v e (t, x) 
Substituting this expression for v into (62) (for a single particle) and we get 



, dV (t ) 
dH 



= q 



e 



e 

dx+^-d t r e (t) x 

c 



B ox dx + e 3 



where 



e 3 = -/ ft 



(x-r f (t)) I x B cx dx 



(67) 



(68) 



V • 



it 



0' 



(x-r*(f)) 



x B fix dx. 



Observe then that according to the charge normalization 



dx = </E cx (t, r e (t)) + ei, where 



(69) 



61 



(E cx (i,x)-E ex (i,r £ (t))) 



0' 



dx, 



and 



B cx dx = B cx (t, r e (i)) + e 2 , where 



(70) 



£2 



(B cx (t,x)-B« (t,r £ (t))) 



dx. 



Hence (67) can be written in the form 



m 



; d V jt) 
d 2 t 



/E cx (t, r e (t)) + 6l + ?-d t r l (t) x [B cx (t, r £ (t)) + e 2 ] + e 3 . (71) 



2.2.1. Dynamics in an external electric field. Here we consider in detail a simpler 
case where external magnetic field is absent, A ox = 0. In this case 



B, 



0, E c 



£2 = 0, e 3 = 0, 



and the field equations take the form 



2m e 



vV 



(72) 
(73) 



As we already mentioned, in this section index £ takes only one value. First, let 
us recall some properties of equations (72), (73). Note that equation (72) does not 
involve the potential ip and can be solved independently from equation (73). If 



14 



ANATOLI BABIN AND ALEXANDER FIGOTIN 



the nonlinearity g{ijj) = G' a (\ip\ 2 ^ ^ is a continuously differentiable function of ijj 
with uniformly bounded derivative, and <p ox is a smooth bounded function, then 
the local existence and uniqueness of a solution of (72) with a prescribed initial 
data in H 1 (IR 3 ) which belongs to and is bounded in H 1 (R 3 ) is well known, see [15, 
Remark 4.4.8]. This solution is defined for alH > if (110) holds for G e since the 
energy 



2m e 



2 



G l 



0* 



2 



dx (74) 



is bounded from below on the invariant set £ H 1 , \\ijj\\ = l}. The existence and 
uniqueness of solutions of the initial value problem for (72) with ip ex = and the 
logarithmic nonlinearity as in (48) were proven in [16], see also [15]. 

Now we will show that, under a localization assumption, the centers r £ (t) of 
wave function if) defined by (58) converge to trajectories governed Newton's law 
of motion for point charges. To eliminate the time derivatives in equation (67) we 
integrate it twice in time obtaining 

r e (t) =4 + r* t+ / ^Ee, (t, / (f x )) dtidta + e 4 , (75) 
Jo Jo m 



where 



1 



m 



o Jo 



£4 = — j / / ei dtidt 2 (76) 



/ V / (E cx (t,x)-E cx (t,r £ (t))) 

JO JR 3 



dxdtidt2- 



1 

m e 

Let us denote z (t) = v l (t) — Tq — fgt and rewrite (75) in the form 

z = F(z)+e 4 , (77) 

F{z)= / -^E ex (t,/(0)+r^(0)i + z) dtid* a . 
Jo Jo m 

We consider now equation (75) on a small time interval [0,T] with solutions in the 
space of continuous functions C ([0,T]). Since the first order spatial derivatives of 
E ox are bounded uniformly, the operator F satisfies a global Lipschitz estimate 

\\F ( Zl ) - F (z 2 )|| c([0>T]) < CT 2 || Zl - z 2 || c([0 T]) . (78) 

In addition to that we have 

\\F(z)\\ cmT]) <C,T 2 if ||z|| c([0T]) <2|r £ (0)|. 

Let us take T > so small that 

CT 2 < 1/4, C{F 2 < 1/4, |r £ (0)T| < \r e (0)| , (79) 

and let Zq be a solution of the equation zq = F (zq) satisfying ||zo|| C ([ t]) — 
2 | (0) | . Such a solution exists and is unique by the contraction principle. Note 
then that the corresponding r (t) = Zq (t) + Tq + r^t is a solution of (75) with e 4 = 
and consequently is a solution of the equation 



p 

m 



^=A x (r), r(0)=ri^m=ii (80) 

which coincides with Newton's equation (66) for a single point charge in an external 
electric field. To estimate a difference between the Newtonian trajectory r (t) and 
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the charge position r e (t) defined by (58) we subtract from (77) equation for Zq and 
obtain 

z-z Q = F{z)-F{z Q ) + e 4 . 
Using (78) and smallness of T we obtain the following inequality 

ll r£ ~ r llc([0,T]) = H Z ~~ Z o!lc([0,T]) - 2 ll e 4|lc([0,T]) ■ ( 81 ) 

The above argument proves the following theorem. 

Theorem 2.1. Suppose that A cx (t, x) = and ip e ^ (t, x) are two times continuously 
differentiable functions, and that the first spatial derivatives of E cx are bounded 
uniformly in (t,x). Suppose also that: (i) there is a family of solutions tp = tp a of 
the NLS (72) with G = G a which depend on the size parameter a > 0; (ii) r £ (t) as 
defined by (58) is a continuous function oft £ [0,T], where T is a sufficiently small 
fixed number which satisfies (79); (Hi) (75) holds where e 4 satisfies the following 
relation 



lim I 

a->0 



\C([0,T]) 



= 0, where \\e 4 \ 



C([0,T]) 



(82) 



sup 

te[o,T] 



^ / (E ox (t,x) 
m ,/ R 3 



E ex (tr f (*))) 



Then ifr(t) is a solution of Newton's equation (80) we have 



sup 

te[o,T] 



1/ (t) - r (t)\ ->0asa 



0. 



Note that the function (E ex (t, x) — E ex (t, r e (t))) = at x = r e and it is a 
bounded differentiable function, implying that condition (82) is a weaker version of 



convergence 



(t, x) -4 8 (x - r e (t)) in C (R 3 ) uniformly on [0, T). Since 

for every a satisfy normalization condition (53), the convergence to 5 (x — r 
means a weak form of localization around r l (t). In Theorem 2.3 and Remark 3 we 
give an example where condition (82) holds. 



(*)) 



2.2.2. Dynamics in a general external electromagnetic field. Now let us consider the 
case when the external magnetic field described by its vector potential A ex does 
not vanish. For localized solutions of equation (30) where N = 1 and £ takes one 
value we want to derive the convergence of r l (t) defined by (58) to a solution of 
Newton's equation of motion 



d 2r „ / x 1 dr 
m— - = gE ox (£, r H rr X B r 



(t,r). 



(83) 



&H 1 cd< 

Considering equation (71) as a perturbation of equation (83) let us estimate the 
difference between their solutions. First we derive from equation (71) estimates for 
(d t r l (t)) similar to usual energy estimates. Multiplying (71) by d t r e (t) we obtain 

m e dy (t) 8 t r e (t) = g f E cx (i, r e (t)) d t r e (t) + (e, + e 3 ) d t v l (t) . 
Integration of the above equation with the respect to t yields 



m 

T 



(sy it)) (o))' 

g £ E cx (t, r £ (<!)) 9 t r f (h) + (e a + e 3 ) «9 t r £ (t x ) dti 

Jo 



(84) 
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Hence 



(d t r e (t)) 2 < (ftr*(0)) +2M 



(ftr £ (h)) dh 



1/2 



where 



M 



< (ft/(0))+Af 2 + / (ftr*(ti)) 2 dt 



dh. 



/ ^( t y(h)) + ^^ 



Consequently, 



(85) 

(d t r e (t)Y < Ah for f G [0,T], (86) 

where constant Mi depends on constants M and T. Let r (f ) be a solution of (83) 
with initial data 

r(0) = r f (0), f£(0) = ^(0). (87) 



(If 



dt 



In a similar way we obtain inequality (86) for (ftr (t)) with a constant Mi = tl/m 
which does not depend on ei + 63. To estimate the difference r e — r we subtract 
(83) from (71), multiply the difference by ft (r £ — r) and integrate the result with 
respect to f obtaining 



I)) 
T 



(ft (r e -r)) 2 = f q* (E cx (i ls r«) - E cx (t, r)) ft (r* - r) df : 

/ ftr x (B cx (t, r £ (t)) - B ox (f, r (f))) ft (r f - r) dh 
c Jo 

+ j ^ft r ^ ( t ) x e 2 + ei + e 3 ) ft (r f - r) dti. 
Therefore, using inequality (86) for (ftr (f)) 2 and (ftr^ (f)) 2 , we obtain 



(d t (r t -r)y<(2C^ 1 + \2CM w 



|ft(r*-r)| dti 



(89) 



+2 J (^ Ml|e2| + |ei+C31 ) l 5 *( r '" r )l dtl 
Using an elementary inequality 



ft (r l - r) dti 



< f 1 / 2 ^ 1 / 2 (t) . 



where 



Z(t)={J Q (ft(r £ -r)) 2 df! 
we obtain from (89) for t G [0,T] 



(90) 
(91) 



(ft(/-r))^< 2C^ + ^V2CMi 
m m c 



+2 





1/2 




dti 







^^(tOlft (r f -r)| dti (92) 

1/2 



<2 1 /2 CT i_ Z (i) + 2Z 1/2 (t) e 5 , 
m 
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where 

\ r T ( q l V 11/2 

Rewriting inequality (92) in the form 

d t Z < 2 1 ' 2 Ct^- 1 Z + 2Z 1 / 2 e 5 , 
mr 

we obtain the following inequality 

z i/2 < ?* , & c,t _ q fortG [o^r], 

with a constant Ci which does not depend on e 5 . Hence Z (t) — > as e 5 — > 0, and 
using inequalities (90) and (92) we conclude that \r e — r| — s- and (pt (r e — r)) 2 — s- 
as £5 -> 0. 

The above argument proves the following theorem. 

Theorem 2.2. Suppose that A ex (i, x) and y ex (t,x) are iwo times continuously 
differentiable functions, and that the first spatial and time derivatives of E cx and 
B cx are bounded uniformly in (t, x) /or i G [0, T] , x G K 3 . Suppose also that (i) there 
is a family of solutions tf> = ip a of the NLS equations (30) with N = 1, G = G a 
which depend on the size parameter a > 0; (ii) r (i) and dtr (t) as defined by (58) 
are continuous functions of t G [0,T]; (Hi) equality (88) holds, and M defined by 
(85) is bounded uniformly in a, < a < ao; (iv) the following limit relation holds 

1 \ 2 

-VMi|e 2 | + |ei + e 3 | ) *->0asa^0, (93) 

77TC / 

where e\,ei, £3 are defined by (69), (70), (68), M\ is a constant which depends on 
M. 

Then if r (i) satisfies Newton's equation (83) we have 

sup |r £ (i) — r (£) | — ^ 0, sup \d t r* (t) — e?tr(*)| ->• as a ->• 0. 

te[o,T] *e[o,T] 

Remark 1. Note that the derivation of Newton's law of motion for charge centers 
in Theorems 2.1 and 2.2 does not depend explicitly on the particular nonlinearity G 
and the value of % > 0. Analysis of the hydrogen atom model in Section 5.1 shows 
that x = h is a natural choice. 

Remark 2. Theorems similar to Theorems 2.1 and 2.2 for a system of interacting 
charges can be analogously formulated and proven. Among additional conditions 
one has to assume the condition of non-collision, that is the trajectories of Newton's 
system (66) have to be separated on the time interval [0, T\. r e (t) — r (t) > d > 

if £ =/= £' with a constant d which does not depend on a. Note that one cannot assume 
that the electrostatic potentials tp have uniformly bounded gradients uniformly in 
a since they converge to the Coulomb potentials. But looking at (82) where E ox 
is replaced by — V<y2 ox — V</S_^ we observe that uniform boundedness of Vip e is 
required only in a vicinity of r e (t) with a fixed radius d/2, whereas outside the 



vicinity a natural assumption is convergence to zero of the integral of 



2 
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2.2.3. Exact wave-corpuscle solutions: accelerating solitons. Here we consider the 
field equations (30), (31) for a single charge, omitting the index £, and present a 
family of their exact solutions in the form of accelerating solitons (wave-corpuscles) . 
Such solutions provide an example for which conditions of Theorem 2.1 hold. We 
assume a purely electric external EM field, i.e. when A cx = 0, E cx (i, x) = 
— V<p ex (i,x). For the purely electric external field the field equations (30), (31) 
take the form (72). We define the wave-corpuscle ip, ip by the following formula: 

ip (t, x) = e iS/x V>, S = mv-(x-r) + s p (t) , (94) 
V> = $(|x-r|), <p = <p(\x-r\), r = r(t). 

In the above formula 4> is the form factor satisfying (35), tp is a radial function 
determined from (31). 

We refer to the function r (t) as wave-corpuscle center. Since ip is center-symmetric, 
this definition agrees with more general definition (58). Note that in a simpler case 
when the external fields tp cx and A ox vanish, a simpler solution of (72) and (83) 
is provided by (94) with r (t) = ro + vt with a constant velocity v. In this case 
the wave-corpuscle solution (94) of the field equations (72)-(73) can be obtained 
from the rest solution ip by certain Galilean-gauge transformations. Solutions 
of a similar form are known in the theory of nonlinear Schrodinger equations, see 
[3G] and references therein. For the particular case of the logarithmic nonlinearity 
solutions of the form (94) were found in [12] in the form of accelerating gaussons. 

Theorem 2.3. Suppose that </J ex (t,x) is a continuous function which is linear with 
respect to x. Then ip defined by (94) is an exact solution to (72), provided that r (t) 
is determined from the equation 



d 2 r ft) 

= <?E CX (t, r) , (95) 



and v (t) , s p (t) are determined by formulas 

v = J, s p = J - qp cx (t, r (i))) dt'. (96) 



Proof. If the potential tp ex (t, x) is linear in x then for any given trajectory r (t) we 
can write 

</?cx (*, x) = ^ ,cx (*) + </4cx (*) ' (x - r (t)) , (97) 

where 

^o,ox (*) = Vx^ cx (r (*) , t) , ip cx (t) = ip cx {t, r (t)) . 
Observe that the representation (94) implies 



dtib = cxp ( i— 
X 



— (v • (x — r) — v ■ r) H - 

X X 



i>-r-Vi>\, (98) 



and by Leibnitz formula we have 
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Substituting the expression (94) into the field equations (72) we obtain the following 
equation for functions v, r, s p : 

[-mv ■ (x — r) — v ■ f — s p ] $ — ixr • (100) 

2 2 

-?v 2 ^ + i X v • V<A + i^V 2 i' - QV e J - = 0. 

2 2m 2m 

Using the charge equilibrium equation (22) we eliminate the nonlinearity G and 
the Laplacian V in the above equation (100) and obtain the following equation 
equivalent to it: 

- jm [v • (x — r) — v ■ r] + y v 2 + s p + ?^ cx | 'tp — i% (r — v) V?/> = 0. (101) 

Now, we equate to zero the coefficients before Vip and tp in that equation, resulting 
in two equations: 

v = f , m [v ■ (x — r) — v • f] + —w 2 + s p + qip eK = 0, (102) 

where, in view of the representation (97), the second equation in (102) can be recast 
as 

2 

to [v • (x - r) - v • f] + s p + ~y~ + q [^ 0jOX + </4cx ' ( x - r )] = °- ( 103 ) 

Then we equate to zero the coefficient before (x — r) and the remaining coefficient 
and obtain the following pair of equations: 

mv 

mv = -qip Qtex (t) , s p - tov • f + + q^o, C x (*) = °- ( 104 ) 

Based on the first equation (102) and the equations (104) we conclude that the 
wave-corpuscle defined by the formula (94) and equations (95), (96) is indeed an 
exact solution to the field equations (72). □ 

Remark 3. The above construction does not depend on the nonlinearity G = G' a 
as long as (22) is satisfied. It is uniform with respect to a > 0, and the dependence 
on a in (94) is only through %p(\x — r|) = a^ 3 ^ 2 ip 1 (a -1 |x — r|) . Obviously, if 
tp (t,x) is defined by (94) then \ip (t,x)\ 2 — >• S (x — r) as a —> and condition (82) 
is fulfilled. 

Remark 4. The form (94) of exact solutions constructed in Theorem 2.3 is the 
same as the WKB ansatz in the quasi-classical approach. The trajectories of the 
charges centers coincide with trajectories that can be found by applying well-known 
quasiclassical asymptotics to solutions of (30) if one neglects the nonlinearity. Note 
though that there are two important effects of the nonlinearity not presented in 
the standard quasiclassical approach. First of all, due to the nonlinearity the 
charge preserves its shape in the course of evolution whereas in the linear model 
any wavepacket disperses over time. Second of all, the quasiclassical asymptotic 
expansions produce infinite asymptotic series which provide for a formal solution, 
whereas the properly introduced nonlinearity as in (22), (36) allows one to obtain 
an exact solution. For a treatment of mathematical aspects of the approach to 
nonlinear wave mechanics based on the WKB asymptotic expansions we refer the 
reader to [24] and references therein. 
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Similarly one can consider a single charge in an external EM field which can have 
nonzero magnetic component. In this case the wave-corpuscle is again defined by 
relations (94) but equations (95) are replaced by (83), see [5] for details. 



3. MULTIHARMONIC SOLUTIONS FOR A SYSTEM OF MANY CHARGES 

In this section we consider a regime of close interaction which models a system 
of bound charges. This regime differs significantly from the regime of remote inter- 
action considered in Section 2. More specifically, we seek special form solutions to 
the field equations (30), (31) for which (i) the potentials ip are time independent; 
(ii) every wave functions ip depends on time harmonically through gauge factor 
e -iutt w jth possibly different values of frequencies u>i for different I. We refer to 
such solutions as multiharmonic. Note that multiharmonic solutions naturally arise 
in the analysis of nonlinear dispersive equations and systems, see [2]. 

Let us consider the field equations (30) and (31) without external fields, namely 



ixdtip 



X 

2m e 



vV 



with ip^ e defined by (23) 



X 

2m l 



1 



<P, if = -q 

The total energy of the system is given by the formula 



l. 



(105) 



(106) 



(107) 



2 



</■' 



dx 



X 

2m e 



dx, 



where ip^ ( is defined in (106), ip £ are given by (32). The energy £ is a conserved 
quantity for the system (105), (106), it can be derived in a standard way from the 
Lagrangian (26) taking in account equations (105), (106). Collecting terms in the 
expression for £ which explicitly involve ip e with a given I we introduce I— th charge 
energy in the system's field by the formula 



-<,)( = 



2m e 



dx. 



(108) 



Here and below integrals with respect to spatial variables are over K 3 . We call so 
defined the £ -th particle energy in system field. Note that for a single particle 
E ( = £( = £, in the general case N > 2 the total energy £ docs not equal the 
sum of E ^, see for example (144) for the case of two interacting charges; in fact the 
difference between the sum of Eq£ and the total energy £ coincides with the total 



energy of EM fields. In accordance with (53) we assume that ip e £ c where 



V- G H 1 



H 2 dx = 1 



(109) 



Theorem 3.1. Let Gg(s) be a function of class C' 1 which is subcritical, namely 
there exists 5 > and C such that 



G, 



dx 



<C + C U\\ 2 H i S for allip€E, £ = 1, .., N. 



(110) 
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Then the functional £ is bounded from below on the set S and boundedness of £ 



implies boundedness of 

' H i- 

to boundedness of all Eo£, £= 1, N. 
Proof. Note that according to (32) 



I = 1, N. Boundedness of £ on S is equivalent 



cp dx 



P (x) 



1 



r (y) 



dxdy. 



|x-y| 

Splitting the domain of integration into |x — y| > 9 and |x — y| < 6 with arbitrary 
small 8 and using the Sobolev imbedding theorem we obtain the following inequality 

2 1 „, 



r (x) 



|x-y| 



■W (y)|^dxdy < 



(111) 



\H 1 



IH 1 



where e can be chosen arbitrary small. From this inequality and (110) we obtain 
the statement of the theorem. □ 

Let us consider now multiharmonic solutions to the system (105), (106), that is 
solutions of the form 



i>t (x) , <P l {t, x) = ip t (x) 



(112) 



Note that we use lower index for ip e (x) to differ it from ip £ (t, x). Then the functions 
ip e (x) satisfy the following nonlinear eigenvalue problem 



X^e^e 



X 

2m t 



V 2 i> £ 



2m e 



-G 



;(km 2 )^ 



0, 



(113) 



where <p^ e is defined by (23), 1 = 1, ...,N. 

3.1. Planck-Einstein frequency-energy relation and the logarithmic non- 
linearity. The nonlinear eigenvalue problem (113), (106) may have many solutions, 
and every solution {V^}f = i determines a set of frequencies and energies 

{Ecwj^Li- Recall that the fundamental Planck-Einstein frequency- energy relation 
reads 

E = huj. (114) 
We would like to find if possible a nonlinearity having the following property: for 
any two solutions {Vvl^Li; {ip'i} »_j with the corresponding frequencies {u!e}f =1 , 
{^>'i\^ =1 and energies {Eo^}^ 1; {Eo^}» =1 the following Planck- Einstein frequency- 
energy relation holds: 

X(^-^) = E ^-E^, i = l,...,N. (115) 

We want the above relation (115) to hold for any regular multiharmonic solution of 
a system of the form (113), (106). Let us consider first the frequencies u>e which can 
be determined from the equations as follows. Multiplying the £-th equation (113) 
by ip\ , integrating the result with respect to the space variable and using the charge 
normalization condition we obtain 



X 

2m t 



IW/dy + gV^r 



dy 



X 

2m e 



G 



i(hk| 3 )hM 2 dy. (116) 



Comparing the above with (108) we see that 

,2 



X^£ 



X 

2m e 



G 



i(l^| 2 )l^| 2 -G/(l^| a ) 



dy. 



(117) 
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For two solutions{^^}^ =1 , we obtain the following relations linking their 

frequencies and energies to the nonhncaritics: 

x(w*-w£)-(EM-Ek)= (118) 

X 2 i si (\ ,/|2\ w (\ ,/|2\ I ,/i2 



Vv| 2 ) - G'i \ |^| 2 dy- 

^(|^| 2 )-GJ(|^| 2 )|^| 2 dy. 



To find a sufficient condition for (115) to hold, we take into account the charge 
normalization condition \\ip(\\ = 1 and observe that it is sufficient to impose the 
following condition 



i (|V>/) - G' e (|V>/) |V>/] dy = const (119) 



for every \tp e \ 2 with = 1; where the constant may depend on £. This condition 
is fulfilled if the following differential equation with independent variable s = 
holds: 

s-^-d (s) - G t (s) = K G s (120) 
as 

with Kq that depends only on I. Then (117) and the normalization condition imply 

2 

G'e (\M 2 ) m 2 ~ Gt (|^| 2 )] dy = ~^jK G (121) 



x 2 

\uj t - E i> = -— 
2m l 



and hence (115) is fulfilled. 

Solving the differential equation (120) we obtain the following explicit formula 
for the logarithmic nonlinearity 

Gt (s) = K G s\ns + Cs. 

Comparing with (50) and (49) we find that if Kg < this is exactly the nonlinearity 
(50) which corresponds to the Gaussian factor. According to (49) the above formula 
takes the form 

Ge (s) = Ge a (s) = o s In s H tts [In _ ,„ — 2 — 3 In a 

y ' ' W (a*f (a0 V ^ 3/2 

where a e is the size parameter for £-th charge. If Kq = then Ge f|Vv| 2 ) ^ s 
quadratic and equation (113) turns into the linear Schrodingcr equation for which 
fulfillment of the Planck-Einstein relation is a well-known fundamental property. It 
is absolutely remarkable that the logarithmic nonlinearity which is singled out by 
the fulfillment of the Planck-Einstein relation has a second crucial property, namely 
it allows a localized soliton solution, namely the Gaussian one in (47). Note that 
above argument can be repeated literally if equation (113) involves external time- 
independent electric field as in (30), in this case (p^ e in the definition (108) of Eoe 
has to be replaced by ip^ e + ip ox . 

Let us study in more detail the situation with the logarithmic nonlinearity singled 
out by the Planck-Einstein relation. Formula (121) for the Gaussian wave function 
takes the form 

x 2 



NEOCLASSICAL THEORY OF ELECTRIC CHARGES 



2:i 



Theorem 3.2. Let every Gg(s), £ = l,...,N, be the logarithmic function with the 
derivative defined by (50) with a = a . Then the functional £ defined by (107) is 
bounded from below on the set S N . The boundedness of £ or boundedness of all 
{Eofl^Li implies the boundedness of and \\ipi\\ X j where the space X is 

defined by (137), for every component of {Vvlfci- 

Proof. The proof of the boundedness from below is similar to the proof of Theorem 
3.1: we use (111), and instead of (110) we use (128). From the boundedness of £ 
or {Eo£j-£ =1 and (128) we infer the boundedness in H 1 , and from the boundedness 



of \G f { \ib,f) dx we infer the boundedness of 



□ 

x 



Theorem 3.3. Let Gi(s), £ = 1, ...,N be the logarithmic functions with derivatives 
defined by (50) with a = a £ . Assume that we have a set of solutions {iplJeLi G > 
a £ £ of the nonlinear eigenvalue problem (113), (106) with the corresponding 
frequencies {o0g}^ =1 and finite energies {Eq^}^. Then for any two solutions tpl = 
ip e and ip" 1 — with a, o\ G £ the Planck- Einstein relation (115) is fulfilled. 

Proof. According to Theorem 3.2 the boundedness of Eq £ for every given a implies 
that functions il>\ are bounded in H 1 and in X. Let us consider the left-hand side 
of (113) and denote it by F (?/v)- Similarly to (111) we obtain that 



tpefVe' dx 



<CUA\m\\ 



iff 1 ' 



where C does not depend on /. Hence multiplication by ip e , is a bounded operator 
from H 1 to H^ 1 . Using this fact to obtain continuous dependence of the term ip^ip^ 
on ip t and using Lemma 9.3.3 in [15] for remaining terms in F(i[i e ), we conclude 
that the functional F (ip#) depends continuously in space W* defined in (138) on 
V-V € H 1 n X and it is bounded in W*. Since 

F (ip t ) Vv^x = Xw* [ IVvl 2 dx - E i+ * — - / \ip e \ 2 dx 
J 2 (cr) rrr J 

(117) holds for smooth rapidly decaying functions ip e , we can use in a standard 
way that smooth functions with compact support are dense in H 1 n X and obtain 
(117) written in the form 

I/ 2 r 

X / I ,<T|2 



x^r-E& = - - ; / i^rdx 

2 (a £ ) m l J 

for functions tp" from H 1 PI X. The right-hand side does not depend on a thanks to 
the normalization condition (53), and we obtain (115). □ 

We would like to stress remarkable universality of the Planck- Einstein relation 
(115) proven above for multiharmonic solutions and the fact that it holds for every 
individual particle in a system of many interacting particles. Namely, the Planck- 
Einstein relation holds with the same coefficient \ f° r an arbitrary pair of solutions 
of (113), (106) and for an arbitrary component of those solutions, and the validity 
of this relation does not depend on the number N which equals the number of 
interacting particles. Observe also that the derivation of (121) and (122) is based 
only on properties of the nonlinear terms as in (118) and can be extended to more 
general systems than (105), (106). In particular, the systems may involve linear 
terms with potentials which explicitly depend on x as in (172). 
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The connection between the logarithmic nonlinearity and the Planck-Einstein 
frequency-energy relation was discovered in a different setting by Bialynicki-Birula 
and Mycielski in [12]. Note though that in [12] a system of N particles is described 
as in the quantum mechanics by a single wave function ip over 3-/V-dimensional 
configuration space, whereas in our model every of N particles has its own wave 
function ip e depending on 3 spatial variables. This is why our approach naturally 
leads to the study of multi-harmonic solutions for interacting particles for which 
every individual wave function ip e can be associated with its individual frequency u)£. 
Another signficant difference between our approach and the nonlinear mechanics in 
[12] is the way the nonlinearity enters the Lagrangian and the field equations for 
N > 2 particles. Namely, in our approach every individual particle has its own 
nonlinearity whereas in [12] there is a single nonlinearity for the entire system of N 
particles. 

3.2. Gaussian shape as a global minimum of energy. As one can see from 
the previous section the logarithmic nonlinearity defined by (48) is singled out by 
the requirement that the Planck-Einstein frequency-energy relation holds exactly. 
Let us look further into the properties of the logarithmic nonlinearity. 

We already observed that the equilibrium equation with the logarithmic nonlin- 
earity admits the Gaussian function as an exact solution. This solution is a critical 
point of the energy, and we show below that the Gaussian function provides a unique 
(modulo translations) minimum of the energy functional on functions subjected to 
the charge normalization condition. This fact was already found in [12]. 

Let us consider the energy £ (ip) = Eq (ip) restricted to one particle with ip ex = 
and G e defined for a = 1 by (49), namely 

£ W = ^I / {|V^ 2 -|^| 2 ln|^| 2 + |^| 2 (ln^-2)|dx. (123) 

Of course, the case a > can be considered similarly. We consider this functional on 
a set 5 defined by (109). In the following theorem we show based on the logarithmic 
Sobolev inequality that the Gaussian function ip g = C g e~l x / 2 provides the global 
minimum of the energy, and that all the global minima belong to the set 

n = {ip : ip = e ld ip g (• - r) , r g R 3 , 6 G R} (124) 

obtained from ip by spatial translations and gauge factor multiplication. Another 
proof of this statement is given in [14] . 

Theorem 3.4. Let ip <E 3. Then, if a = 1 



£ (vo > £ = 




(125) 



and the equality holds only for ip G Q. 

2 

Proof. For simplicity we set here &j = 1. The value of the energy £ (ip) on ip g can 
be found explicitly: 

S{%)=\ (126) 
Note that the Euler equation for constrained critical points of £ (ip) has the form 

+ V 2 iP = (- In (|V;| 2 ) + ( In -ij - 3^ (127) 
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and ip g satisfies this equation. Now we show that tp g provides the global minimum 
of £ (ip) defined by (123) for real-valued ip. We will use the following Euclidean 
logarithmic Sobolcv inequality 



ib In ib dx < — In 

Y Y - 2 



4- / IW| 2 dx 
nae J R d 



(128) 



for functions from H 1 (K d ) where equality holds only for Gaussian functions (see 
[11]). Therefore, 

£{i>)> (129) 



|VV| 2 + H 2 (In 



-3/2 



- 2 



dx In 

2 



-H- f |V^| 2 dx 



-i{/|V*|'dr-l-|h[|^|V*|«fc 

To find J Rd \S7ip\ 2 dx we consider the Euler equation (127) for the minimizer ib. 
From Pohozhacv identity applied to (127) we obtain: 



1 f 1 



-G 



(|*| 2 ) dx = AM 



WW dx, 



(130) 



where 



G 



(Hi 



■ln|^r + |^r In 



r 3/2 



Multiplying equation (127) by tp and integrating the result we get 

A||V,|| 2 -£ W = -\ J kf dx = -i 
Adding (130) and (131) we obtain 



(131) 



X 
'2m' 



J | V? M 2 dx=-i-i J |W| 2 dx, 



and hence 



J \Vtjj\ 2 dx = 



Substitution of the above expression into inequality (129) implies £ (tp) > \ where 
the equality holds only for Gaussian functions. Since tp is a minimizer, the inequality 
holds for any real tp € H. Let us consider now a complex minimizer ip = u + iv such 
that £ (tp) < |. Notice that the following elementary inequality holds 

V(u 2 + i> 2 ) 1/2 | 2 < \\7u\ 2 + \\7v\ 2 
if vr -h v 7^ and the equality is possible only if v\7u = uVv. The above inequality 

1 /2 

implies in particular that (u 2 + v 2 ) £ iJ 1 . Applying once more the Euclidean 
logarithmic Sobolev inequality for real functions we get 



,V2 



2\!/2 



For such a minimizer (u 2 + V 2 ) must coincide with a Gaussian function and 
equality vVu = u\7v must hold a.e.. Since 



Vu Vw „, 
V In w = = = V In v 

u v 
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for non-zero u and v, we conclude that In it — lni> is a constant and the complex 
minimizer ip is obtained from a Gaussian function via multiplication by a gauge 
factored □ 

The following theorem directly follows from results of Cazenave and Lions [17]. 

Theorem 3.5. Let fl e be e -neighborhood in H 1 norm of O 



e = WeE: inf||^ -v\\ H1 <e}, (132) 

|_ v&U J 

and Ms be a sublevel set 

M s = {i/>eEi£{il>)<e(l> g )+5}. (133) 
Then for any e > there exists S > suc/i that Ms C il e . 

Proof. Assume the contrary, namely that there exists e > and a sequence ipj G 3 
such that £ (tpj) — > 5 (ip g ) and inf„ e n ||^ — v || x > e. According to Theorem II. 1 
and Remark II. 3 in [17] a subsequence ipj (• — r 3 -) is relatively compact in H . Hence 
there is a subsequence that converges in H 1 to a global minimizer which belongs to 
fi. This contradicts the assumption e > 0. □ 

Let us consider now one-particle equation (105) with logarithmic nonlinearity in 
the absence of external fields (a free particle) 

2 2 

[ xdt4> = -^ V V + ^ (- ln (H 2 - 3 ) ^ ( 134 ) 

The existence and the uniqueness of solutions to the initial value problem for this 
equation is proven in [16]. Since the logarithm has a singularity at ip = 0, to describe 
classes of functions for which the problem is well-posed we need to introduce special 
spaces. Following Cazenave [15] we introduce functions 

A (\4>\) = - IV>| 2 ln|V>| 2 if < \i>\ < c~ 3 , (135) 
A(\ip\) = 3 |-0| 2 + 4e- 3 |-0| — e~ 3 if |V>| > c" 3 , 

B(\iP\) = \iP\ 2 ln\iP\ 2 + A(\iP\). (136) 
So defined A is a convex C 1 function which is C 2 for |^| ^ 0. The function B {\ip\) 
vanishes for small and tpB (\ip\) / \ip\ is of class C 1 and satisfies uniform Lipschitz 
condition. Let A* be the convex conjugate function of A, it is also convex of class 
C 1 and positive except at zero. We introduce Banach spaces 

X = {V e L 1 (R 3 ) : A(|V>|) G i 1 (R 3 )} , (137) 
X' = {V € L 1 (R 3 ) : A* (|V|) G i 1 (R 3 )} . 

Let 

W = i? 1 (R 3 ) n X, W* = H 1 (R 3 ) + X'. (138) 

Note that Ms is bounded in W. By Lemma 9.3.2 from [15] the function ip — > ip ln l^ 2 
is continuous from X to X' and is bounded on bounded sets. 

According to Cazenave [14], [15], if ip G W there exists unique solution ip G 
C (R, W) n C 1 (R, W*) of (134) with initial data ip (0) = ip . Properties of solutions 
of (134) are described in [14], [15], in particular S n W and Ms are invariant sets. 
According to [17] the solution ip is orbitally stable, namely invariance of Ms and 
Theorem 3.5 imply that if ip (0, •) G Ms then ip (0,t) G f2 e for t > 0. Consequently, 
small initial perturbations of the Gaussian shape do not cause its large perturbations 
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in the course of time evolution given by (134) but may cause considerable spatial 
shifts of the entire wave function. Such a time evolution of charges in electric field 
is described in Theorem 2.1. Note that if the potential ip eK (t,x) is linear as in 
Theorem 2.3, one can rewrite equation (72) in a moving frame with origin r(t) 
and corresponding phase correction for tfj as in (94) (see [5] for details) and obtain 
an equivalent equation of the form (72) with ip ex (t, x) = 0, namely (134). Wave- 
corpuscle solution (94) after the change of variables turns into ijj g . Therefore, using 
Theorem 3.5 we conclude that the wave-corpuscle solutions of (72) constructed in 
Theorem 2.3 are orbitally stable. 

4. TWO PARTICLE HYDROGEN-LIKE SYSTEM 

In this section we consider a particular case of multiharmonic solutions of the 
system (113) with N — 2 which models a bound proton-electron system. For this 
system we discuss general properties and provide a heuristic analysis of the coupling 
between proton and electron. Using this analysis as a motivation, wc then write 
equations and energy functional for one charge (electron) in the Coulomb field of 
proton. We study this problem in more detail along the lines of [9], [10]. 

4.1. Electron-proton system. To model states of a bound proton-electron sys- 
tem, the hydrogen atom, we consider the multiharmonic solutions of (30), (31) 
described by the system (113), (106) with N = 2, where indices £ take two values 
I = 1 and 1 = 1. The charges have opposite values q\ = —q 2 = —q. For brevity we 
introduce notation 

$i = tpjqi, $2 = V2/92, (139) 

ai = -f— , a 2 = . (140) 

g 7711 q 777,2 

The quantity ai turns into the Bohr radius if x coincides with Planck constant H, 

and mi , q are the electron mass and charge respectively. Using the above notation 

we rewrite the two-particle system (113), (106) in the form of the following nonlinear 

eigenvalue problem 



2^ lV)l + ^Lv 2 ^ + $2^ = y g 'i (V'il 2 ) ^1. ( 141 ) 

4^2^ 2 + y V 2 ^ 2 + $1^2 = yG' 2 (|^ 2 | 2 ) ^2, (142) 



V 2 $x = -47r|V>i| 1 V 2 $ 2 = -477^2^ ■ (143) 
The functions ip 1 and i/) 2 are respectively the wave functions for the electron and 
the proton. As always, we look for solutions ip t £ 5 with 3 defined by (109). 
The nonlincarities G^, G' 2 are assumed to be logarithmic as defined by (50) where 
the size parameter a = a 1 is different for electron and proton. One can similarly 
consider other nonlinearities but in this paper we stay with the logarithmic ones. 

In accordance with (107) we introduce the energy functional by the following 
formula 

£(1>i,1> 2 ) (144) 



e J [y iVV'rl 2 + f G X (|^| a ) - 4. |^| 2 (-V 2 )" 1 |^ 2 | 2 



dx 
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where (-V 2 ) 1 |Vv| 2 is defined by (32). Equations (141), (142), (143) can be de- 
rived from the Euler equations for a critical point of £ (il)\,tp 2 ) with the constraint 



(^, V 2 ) G S 2 , if we set in (143) $ 2 = ^ [(-V 2 ) 1 |^ 2 | 2 J , $1 = 4tt ((-V 2 ) x \i/j 2 \ 
the frequencies uii,lu 2 being proportional to the Lagrange multipliers. Using (32) 
one can see that 

J i^ 1 i 2 (-v 2 )- 1 i^ 2 i 2 dx= J iv/t-vTV^dx 

hence the coupling term in (144) is symmetric with respect to ip\,'4>2- Note also that 
the first term in formula (144) coincides with the energy E01 of the first charge in the 
system field as given by (108) with £ = 1, and the second term with the expression 
for the energy £ 2 {ip 2 ) °f a f ree second particle given by (123). Obviously, £ ip 2 ) 
also can be written as a sum of energy E02 of the second charge in the system field 
plus the free energy of the first charge. 

According to Theorem 3.2 the energy £ (ip 1 , tp 2 ) 1S bounded from below. 

Remark 5. Note that system (141)-(143) is similar to the Hartrcc equations studied 
in [27], though it differs from it because of the presence of the nonlinearities Gg. 

4.2. Reduction to one charge in the Coulomb field. In this subsection we 
give a motivation for replacing the problem of determination of critical values of 
the energy functional £ given by (144) for a system of two charges by a simpler 
problem for a single charge similarly to the Born-Oppcnhcimer approximation in 
quantum mechanics. To this end we use two changes of variables in two equations, 
x = aiyi in (141) and x = a2y2 in (142), where at, are defined by (140), namely 

x = a t yi, f = 1,2, (145) 

and rescalc the fields as follows: 

^(x) = ^M ^(x) = -±a**(y*)>*=l>2- (146) 

a e a/ 

Hence (141), (142) takes the form of the following nonlinear hydrogen system 

*i| 2 )*i, (147) 

G' 2 (|* 2 | 2 ) #2, (148) 

V 2 ^ ^^I^l 2 , V 2 2 2 = -4tt|v& 2 | 2 . (149) 
Here we use the same letter G to determine the function of rescaled variables. Note 
that the electron/proton mass ratio ^ ~ is small, therefore the parameter 

b=^ = ^«l (150) 
ai m 2 

is small too. We rewrite (147), (148) as follows: 

xai _ , i„ 2 ,t, , 1, (y 



72\-! |„/. |2 



q 


1 

~ 2^ 


V 2 x *H 


"-02 

a 2 




*i = 


1 

2a~! 


4^2*2 "1 

Q 


1 

" 2a~2~ 


V 2 2 * 2 H 


--01 
ai 


(a 2 > 
— Y2 


*2 = 


1 

~2a~2 



<1 2 



+ + \<t>2 (|) *i = \g\ (l*i| 2 ) *x, (151) 



^w 2 *2 + ^V 2 * 2 + Hx (by) *2 = \g' 2 (|* 2 | 2 ) *2- (152) 

We would like to show that if we are interested in electron frequencies u)\ which 
correspond to lower energy levels of £ (tpi, ?/> 2 ) then, in the case b — > 0, it is reason- 
able to replace system (151), (152) by one equation with the Coulomb potential. 
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The discrete spectrum of the classical linear non-relativistic Schrodinger operator 
for hydrogen atom with the Coulomb potential can be recovered (if multiplicities 
are not taken into account) from radial eigenfunctions. Therefore here we restrict 
ourselves to the radial solutions of (151), (152) and restrict (144) to radial func- 
tions. Consider a radial solution (vf^, ^2) of (151), (152). Equation for the electron 
frequency uj\ given by (116) with I = 1 takes the form 

xai. . _ f ( i,t, ,2\ ,, T , ,2 1 IV7lT , ,2 A , /y\ 1 \ ,, T , ,2 1 ,, T . ,2 



(153) 

with the proton potential <p 2 . Then we would like to replace the potential cf> 2 in 
(153) by the Coulomb potential t^t as an approximation. To show that for small 
b one may expect the resulting perturbation of eigenvalues to be small, we use the 
following representation for the radial potential <j> (r) = An (—V 2 ) \^\ 2 - 



4>{r) = - 



1 - An I (n - r) ri |* (r x )| 2 dri 



(154) 



The term which we expect to be small is 



[\<t>2 (JyiJ - rMl*i| 2 d 



= An I — I*! (r)| 2 r 2 / (n - r/6) rj |* 2 (r^ 2 dndr. 

^0 r Jr/b 

Changing the order of integration and using substitution r/b = r 2 we obtain 

/>oo />oo 

D pro t = (Anf b 2 / |*i(&r 2 )| 2 r 2 (r 1 -r2)r 1 |*2(ri)| 2 d n dr 2 

JO Jr 2 
/•oo />ri 

= {Anfb 2 \^> 2 {r 1 )\ 2 r 1 (6r 2 )| 2 r 2 (n - r 2 ) dr 2 dr x . 



Hence 



l^protl < ^& 2 max|\E'i(r)| 2 /°°|* 2 (ri)| 2 r4dr 1 . (156) 
6 Jo 

The right-hand side of (156) involves the variance of if? 2 and H^iH^no ■ To show 

that |-Dpiot| is small it is sufficient to show that the variance of ^2 and H^iH^oo 

are bounded uniformly for small b. Note that ^1 , vE r 2 in (156) are not arbitrary 

radial functions, but special solutions of (151), (152), (149) which correspond to 

lower critical energy levels of the energy functional £ = E i + E 2- Dependence of 

the coupling between (151) and (152) on b comes through the positive potentials 

\ < t ) 2 i\y) an d b(j) 1 (by) in (151) and (152) respectively, according to (154) these 

potentials are bounded by A uniformly for b > 0. Using the uniform boundedness of 

the potentials and estimate (182) we conclude that energy and H^-norms of \&i and 
\I/2 are bounded uniformly in b. As in the proof of Lemma 5.3 below, we conclude 
that radial "Ji and if? 2 arc bounded in L 5 , they also satisfy (210). From boundedness 

in L 5 we derive that j^^i is bounded in L 7 / 4 and so does G[ fl^il 2 ) ^i- Note 



A"l . 
q 2 

in (151) which correspond to lower energy levels are bounded. Hence, the term 



also that since the energy is bounded from below, the negative values of u> = ^r^\ 



b^\b J l- 2 



1 \ T 1 



y *i--Gi |*iP)*i + o;*i 
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in (151) is bounded in L 7 / 4 , therefore the solution "fi of the elliptic equation (151) 
belongs to Sobolev space W 7 ^ 4 ' 2 and by Sobolev imbedding is bounded in L°°. As 
for ^2, from (210) using the comparison argument as in the proof of Lemma 5.8 we 
derive exponential decay of solutions of the equation (152) for large r uniformly in 
b and conclude that the variance of ^ is uniformly bounded for small b. Hence, for 
such solutions |-D pro t| is small for small b. The smallness of b in (156) (recall that 
b ~ ygTjy for proton-electron system) gives us a motivation to set b = in (152) and 
replace \(f> 2 (±y) by ^ in (151). 

So, we replace the problem of finding frequencies cjj by formula (153) based on 
critical points of £ (ip 1 ,i[i 2 ) which correspond to lower energy levels of £ ip 2 ) 
by the problem of finding critical points, lower critical levels and corresponding 
frequencies wi (which for the logarithmic nonlinearity are expressed in terms of 
the critical levels by (122)) for the following energy functional with the Coulomb 
potential: 



£ 



Ob 



| |V*iI a + id (|*i| 3 ) - I*il a 



dy. (157) 



This functional involves only ^i. In quantum mechanics a somewhat similar re- 
duction (in a completely different setting) from proton-electron system to a single 
equation with the Coulomb potential is made via Born-Oppenheimer approxima- 
tion. Additional motivation for the reduction is given in Remark 10. 



5. Variational problem for a charge in the Coulomb field 

In this section we consider in detail radial critical points of the functional £cb (*i) 
defined by (157) and establish its basic properties. In this and following section we 
look for real solutions ^ — *!/* and all function spaces involve only real functions. 
In what follows we explicitly take into account the dependence of the nonlinearity 
on the size parameter a = a 1 (or on related parameter k), namely 

G' /r (\*\ 2 )=k 2 G' 1 (k- 3 \*\ 2 ), 

where, for consistency with the notation (33) where a = <Xi/k, we use notation 
/ k to denote the dependence of rescaled G on k. Note that the dependence on a\ 
and a = a 1 is factored out after the rescaling, and that G/ K in the above formula 
depends only on their ratio k. We consider the case of small K, namely 

k= — <1, (158) 
a 

when the electron size parameter a is much larger than the Bohr radius a\. To 
treat technical difficulties related to the singularity of logarithm at zero, together 
with the logarithmic nonlinearity given by (50) we use a regularized logarithmic 
nonlinearity defined as follows: 

G% l£ (lV'1 2 ) = W,e ( K ~ 3 l c t) ~ 3k2 > ( 159 ) 

where £ < 0, 



ln +! £ (s) = max (In (s) , f ) for s > 0, ln +i? (0) = £. 



(160) 
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For £ = — oo we set ln+^oo (s) = ln(s). The function ln + ^ (s) with finite £ is 
bounded for bounded s and satisfies Lipschitz condition for s > 0. Obviously, 

- ln + , 5 (s) < - ln+, r (s) for £ > £' > -oo. (161) 

Now we describe the basic properties of G'^ ? and its integral Gy re ,g 

G hi (s)= [ S G' lti (s')ds', G /Kti (s) = - K 5 G U (^ 3 S ) . (162) 
Jo 

Proposition 5.1. T/ie function G/ K ^ with — oo < £ < can &e written in the form 
G /k4 (s) = - K 2 S (ln + , e {n- 3 S /C 2 g ) + 2 + C ff 2 c«-'»+, 5 



where 



e £-W,«« < m in [1,— ] < 1 



(163) 
(164) 



and satisfies the inequalities 

G/k,$ (s) < G /K , e (s) < G/^ (s) = ~k 2 s In (a) - n 2 s (in jC] + 2) , (165) 
for > £_ > £' > -oo, 



Gi, S f|*| 2 ) - Gi,_oc (|Mf) <Gc«/ 4 |*| 



Proof. Elementary computation shows that for s > G 2 e £ 

G u ( S ) = -sin ( S /G 2 ) - a (C^" 1 " + 2) 



and 
Hence 



G u 0) = -C- s - 3s for s < e? C; 



G U ( S ) = -s (ln + , 5 ( S /G 2 ) + 2 + G 2 c«- ln -« s ) 
and (163) holds. Using (161) we derive (165) from (162). Obviously 
G u (s) - G^-oo (s) = -e £ for s > e«G 2 , 

and hence 

G u (|*| 2 ) - Gx.-oo (|*| 2 ) | < e«/ 2 |*| /G 9 for \V\ 2 > c«G; 
For |*| 2 < e«G 2 

Gi,e (|*| 2 ) - Gi,-oo (|*| 2 ) I = |*| 2 I (£ - In |*| 2 /G 2 
< e« /4 G 9 |*| 3/2 In |*| 2 /G 2 < Cie«/ 4 G 9 |*| . 
^Frorn (170) and (171) we obtain (166). 



(166) 

(167) 
(168) 
(169) 



(170) 



(171) 



□ 



The nonlinear eigenvalue problem with the Coulomb potential has the form sim- 
ilar to (151): 



1 



1 1 



w * + _ V 2 * + — vf = -G' /Kti [\^) *, 
where the dimcnsionless spectral parameter ui = uj k ^ is given by the formula 



(172) 



Xai ,i 



(173) 
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Note that equation (172) can be obtained as the Euler equation by variation of the 
energy functional (157) which we write in the form 

= £b (*) + (*) , ( 174 ) 

where £ o is the quadratic energy functional 

£0= f (l\W\ 2 -^\*\ 2 ) dx, (175) 



2 1 1 |x, 
and Qk,c, is the nonlinear functional 

g K , £ (*) = - y / (in+, 5 («- 3 l*l 2 /c 2 ) + 2 + «- 3 i*i 2 ) |*| 2 dx. (176) 

As always we assume the charge normalization constraint 

||*|| 2 = / |*| 2 dx = 1 (177) 



and denote similarly to (109) the set of radial functions (which depend only on |x|) 

by 

2rad = {*eH I ] ad (K 3 ):||*|| 2 = l}. (178) 
Obviously, for $ € H n X where X is defined by (137) 

& (*) = &,-«,(*) = — J (|*| 2 ln|*| 2 ) dx+y (in -L-2 + 31m 

(179) 

The spectral parameter oj is the Lagrange multiplier and it relates with correspond- 
ing critical energy levels of £ K ^ on S fl X by formula (122) which takes the 
form 

= w„, f>B + y. (180) 

Let us introduce for K > 0, £ < 0, \& S C the function 

(*) = (I*! 2 )* = (-K 2 ln+, e (V 3 |*| 2 /C 2 ) - 3k 2 ) *. (181) 



General properties of functional £ K .£ ( , 5) with regularized logarithm are described 
in the following statements. 

Lemma 5.2. The functional £ K £ (\1/) defined by (174), (169) and (176) with —00 < 
£ < on S ra( j /ias the following properties: (i) it is bounded from below on H 
uniformly in £. (ii) f/feS, k<1 and (*) < C i/ien < C" where C 

depends only on C. (Hi) The functional £ K ^ (\f) is of class C 1 with respect to H 1 
norm on S ra( j. (iv) 

G K ,e (*) > Qn,o (*) > -C"« 2 (1 + |in«|) if \\n H i < C". 

Proof. Boundcdncss from below follows from (129), (165) and inequality 

^l^l 2 dx< Co ||*|| ||V*||. (182) 

R3 |X| 

^From boundcdncss of £ K ^ (*), (129) and (165) we obtain boundedness of ||V\&|| 
which on S implies boundcdncss in H , To prove (hi) observe that by (159) 



ff «, f (*) = K 2 (-e-3)*+. i 5 (|vi/| 2 ) 
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where g\^ {^f^ is identically zero for small and it has less then quadratic 

growth for |*| oo. According to Theorem A.VI in [9] then Q K £ (*)+k 2 (£ + 3) ||*|| 
is of class G 1 , and hence Q K £ {^) is of class G 1 too. To obtain (iv) we use (165) 
and observe that Q K) q (\I/) can be estimated in terms of H^H^e which, in turn, can 
be estimated by H^H^i. □ 

Lemma 5.3. Let a sequence of radial functions € S r ad satisfy (183). Then: (i) a 
subsequence of ^ j converges in LP (R 3 ) and almost everywhere to if 2 < p < 6; 
(ii) the functions vfj (x) — > as |x| — > oo uniformly in j. 

Proof. The estimate (183) together with the normalization condition imply the 
uniform boundedness of \£ j in the Sobolev space H 1 (R 3 ) . According to Theorem 
A.I' of [9], the space of radial functions H^ ad (R 3 ) is compactly imbedded into 
LP (R 3 ) if 2 < p < 6. Hence, a subsequence (x) -> *oo (x) strongly in LP (R 3 ) 
and almost everywhere. Statement (ii) follows from Radial Lemma A. II in [9]. □ 

The functional £ K .^ (VP) with regularized logarithm satisfies the Palais-Smale con- 
dition similar to condition (P-S + ) in [10]: 

Theorem 5.4. Let £ > — oo. Let a sequence £ S ra d have the following proper- 
ties: 

£«,f(%)< -P, /?>0, (183) 

and 

F* s = + iv 2 *, + -L^ - lg Kti -> (184) 

strongly in H~ 1 (R 3 ). Then the sequence contains a subsequence which converges 
strongly in H^ ad (R 3 ) . 

Proof. The proof is similar to treatment of radial solutions in [9], [27]. We only 
sketch main steps. Since (182) holds, (183), (165) and (128) imply boundedness of 
tyj in H 1 (R 3 ). Using Lemma 5.3 we conclude that we can choose a subsequence 
such that _> vp^ in L 5 (R 3 ) and L 8 / 3 (R 3 ) , it converges weakly in H' 1 (R 3 ) 
and almost everywhere in R 3 . Multiplying (184) by in L 2 (R 3 ) = H with inner 
product (■, ■) we obtain 

(*i.*j> - \ (V*i,V* 3 -) - \{(G)^ I*/) U-^) (185) 

We rewrite (185) in the form 

Wj - = (1?*^,) + (*,•) - | G /K , ? (|vl/| 2 ) - G' /K ^(|vl/| 2 ) |vl/| 2 dx. 

Using (164) we obtain for \I> £ S rac j 

< j G/^(|*| 2 ) - G^(|*| 2 ) |*| 2 dx < ^. (186) 

Since (F^j^j) — > 0, (183) and (186) imply that we can choose a subsequence tyj 
with u>j — > where < — J3. Multiplying (184) by a smooth test function 
with compact support and passing to the limit we obtain in standard way that 
*oo G H 1 (R 3 ) satisfies (172) with lj = uj oq . 
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Now we want to prove that ^oo G S and convergence in H^ ad (R 3 ) is strong. We 
rewrite (185) in the form 

- \ (V^-,V%) - \ (|*,| 2 ) - K 2 e) (187) 



1 



We have (FV^Vj) 0, 



Since G K .£ (J\&| J — k £ is identically zero for small I'J'I and grows slower than |\I>| 

as 1^1 — ¥ oo, and are converging in L 1 (R 3 ), we can apply the compactness 

Lemma of Strauss (see Theorem A.I in [9]) and obtain 



G 



/«,« (i*,f) -« 2 e) -+ ((g /k , c (Vooi 2 ) -« 2 e) *oo,*c 



Using the Holder inequality we obtain that the quadratic mapping \I/ — ^ jij l^&l 2 is 
continuous from L 5 (R 3 ) n L 8/3 (R 3 ) into L 1 (R 3 ) . Therefore 



We denote 



M = lim inf = 1, Mi = lim inf (V^.V^) . 

j—¥OC j—tOQ 

^From (187) we obtain passing to the limit 

( Woo + K 2 i) M - ^Mi - l((G /KA (|*oo| 2 ) - k 2 *) *co, *co) (188) 

1 

l x l 

Multiplying (172) with cj = by in H wc obtain 



Woo (*oo,*oo) 



i (Woo, Vv&oo) - \ ( (G/^ (|*oo| 2 ) - K 2 C) *oo, * c 



Comparing with (188) we obtain 

(Woo + K 2 C) Mo - ^-hh = (Woo + K 2 ^ (*oo, *oo) - 2 (Woo, Woo) . (189) 

^From weak convergence in ff^d (R 3 ) we infer that M > ||^oo|| 2 , M x > ||Woo|| 2 - 
Since (woo + k 2 £) < 0, (189) is possible only if M = |^oo| 2 , Mi = |Woo| 2 - These 
equalities imply that the weak convergence of the subsequence in H 1 to Xt'oo is 
strong. □ 
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5.1. Nonlinear eigenvalues for a charge in the Coulomb field. In this sec- 
tion we prove that if k is small enough, the nonlinear eigenvalue problem (172) 
with logarithmic nonlincarity has solutions with eigenvalues which are close to the 
eigenvalues of the corresponding linear problem. We look for real-valued solutions. 
As a first step we prove existence of such solutions of the problem with regularized 
logarithm ln +j £. As the second step we pass to the limit as £ — > — oo and obtain 
solutions of the eigenvalue problem with the original logarithmic nonlincarity. 

Consider the linear Schrodinger operator O with the Coulomb potential which 
corresponds to £ (^)j that is 

0* = -^V 2 * - |-|*. (190) 
2 |x| 

Note that O has the following well known negative eigenvalues 

E Q n =cj , n = -^, n=l,2,... (191) 

which coincide with the negative eigenvalues of the operator O ra d obtained by re- 
striction of O to radial functions (without counting their multiplicity) and are equal 
to the negative critical energy levels of £q (VP). Every eigenvalue of the radial prob- 
lem 

is simple and corresponding eigenfunctions (r) given by well-known explicit for- 
mulas decay exponentially as r — > oo, for example the linear ground state is given 
by the formula 

Let us denote by lP n an eigenspace of r ad which corresponds to the eigenvalue 
E® and consists of radial functions, by L~ C H^ ad a finite-dimensional invariant 
subspace of O ra d which corresponds to eigenvalues A < E® and by L+ C H^ ad (M 3 ) 
invariant infinite-dimensional subspace in H^ d (M 3 ) of ra d which is orthogonal in 
L 2 to L~ . We have an orthogonal decomposition 

H Ld ( R3 ) = L n + L n = L n-1 + L n + L n- 

We look for eigenvalues u> = cj k t n < of the nonlinear eigenvalue problem (172) 
which can be considered as a perturbation of the eigenvalues E®. The ground state 
of £ K ,-oo (*) and the corresponding energy level w«,-oo,i was found in [12], it is 
given by the formula 

*«,-oo,i (x) = K 3 / 2 C s c-l KX l 2 / 2 Cc-l x l, = - J - \k 2 InC 2 , (192) 



J g 
where 

c- 2 



r l/2 



^V^e- ^v^erf (l)e- - 



2 2 
iln C 2 cli 



Hence the nonlinear correction for the ground level is of order k 2 . 

Since the Palais-Smalc condition is satisfied according to Theorem 5.4, the ex- 
istence of an infinite sequence of discrete negative energy values of £ K ^ (*) can be 
proven exactly along the lines of [9]. But we want to show, in addition, that the 
eigenvalues of the nonlinear problem with small k are close to the eigenvalues of the 
linear problem (this is not surprising, since all the eigenvalues of the linear problem 
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for radial functions arc simple). Therefore we modify the construction of critical 
points of the functional £ K ^ on S. 

We will use the following deformation lemma which is a minor modification and 
reformulation in our notation of Lemma 5 from [9] (with essentially the same proof). 

Lemma 5.5. Suppose that the functional £ K ^ £ C 1 (S ra d,R) satisfies Palais-Smale 
condition as in Theorem 5.4- Suppose also that E" < and that the segment 
E' < E < E" does not contain critical values of £ Kt £ . Let e > be arbitrary small 
number. Then there exists < e < e and a mapping r\ g C (H ra( j, H ra d) such that 
(i) i] (u) — u for u G S ra d with £ K .^ (u) < E' — e or £ K £ (u) > E" + e. 
(ii) rj is a homeomorphism S ra d — > S ra d and it is odd if £ Kl £ is even. 
(Hi) £ K £ (rj (u)) < £ Kii (u) . 

(iv) r\ ({it g H rad : £ K £ (u) < E + e}) C {u e S rad : £ K £ (it) < E - e} for E g 
[E',E"]. 

To prove existence of critical points of £ K ^ ( , J) on the infinite-dimensional sphere 
Srad with critical values E^ near E%, n = 2,3... when k is small, we slightly 
modify the construction from [9] . As a first step we define finite-dimensional domain 
B = B n C L~ which lies in S ra d near eigenfunction vf^ of O ra d as follows: 

B n = (1 - R 2 ) 1 ' 2 *° +v, ve L-_i, ||v|| = R, 0<R< R ] , (193) 

where Ro < 1/2 is a small fixed number, for simplicity we fix Rq = 1/3. The domain 
B n is diffcomorphic to a n — 1 dimensional ball. The boundary of B n is the n — 2 
dimensional sphere 

dB n = {^ :V=(l-I%) 1/2 * n + v,veL-_ 1 , \\v\\=R }. (194) 

We consider continuous functions T defined on B n with values in S. We consider 
the class {r}„ of continuous mappings T : B n —> S ra d such that T restricted to dB n 
is identical, 

r(#) = * if * g dB n . (195) 
The n-th min-max energy level E K ^ <n is defined as follows: 

E K ^ n = inf max^(r(f)). (196) 

Lemma 5.6. Let £ < 0, let n = 2,3,...., < K < 1. TTiera E K £ n defined by (196) 
lies in the interval 

I n = \E',E"} , with E' = E°- C 3nK 2 (1 + |ln«|) , E" = E° + C ln K 2 , (197) 
where C\ n and C 3n depend only on n. 

Proof. Since we can take a particular function V ("J) = To (^>) = + v for g B n 
we have 

K,t, n < max £ K , e (*) < £ + max^ (*) = £° + max£ K , C (*) ■ (198) 

'I'G-Bn B n B n 

Note that B n lies in a ball in a finite-dimensional subspacc L~ in H 1 , the basis 
of L~ consists of exponentially decaying functions with bounded first derivatives, 
hence B n is compact in H 1 (R 3 ) and in L p (R 3 ) for every p > 1. Therefore 

+ ll*llia + < C 0,n for * G B n . 

Hence, using (165) and power estimates for ^ 2 ln\I' 2 we conclude that 

(*) < a K ,-oc (*) < C« 2 (||*|| 3 L3 + ll*!!^) - K 2 (\nn- 3 /C 2 + 2) ||v|/|| 2 
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and we obtain 

maxg Kt6 (*) < C ln K 2 , (199) 

where C\ n depends only on n. Using Tq (iff) and the estimate (198) of £ K £ (iff) we 
observe that 

E K , $ , n = inf max £„, f (T (¥)) , (200) 

where C 2n = Ci„k 2 + E°. Inequality (199) provides E" in (197). 

Now we estimate E^^^ n from below. Mappings T € {r}„ which satisfy (195) 
have the following topological property: the set T (B n ) must intersect the subspace 
L^-i- Indeed, assume that the contrary were true and for some T there were not 
intersection of T (*) with that is UlL^r (*)|| ^ for all * g B n where 

n~_j is the orthoprojection onto the finite-dimensional subspace L~_ x . Then we 
could consider the following function on the ball ||w|| < Ro in L~_^. 

T l (v) = ^ r) U-_ 1 T( V ), 

where N (v) = n^_i (r (v))\\ . This function is continuous on the ball {||i>| < Ro} 
(since the denominator N (v) does not vanish) and Li (v) — v if ||w|| = f?o, and also 
it has the property \\Ti (v)\\ = Ro if \\v\\ < Ro- A function with such properties 
is called a retraction of the ball {\\v\\ < Ro} in the (n — l)-dimcnsional space L~_ x 
onto its boundary{||i;|| = Ro}. It is a well-known fact from the algebraic topology 
that such a retraction cannot exist. Therefore for every T we have II^'_ i r ('I'o)ll = 
for some^o € B. Therefore there exists V^o € B n such that 

r (*(,) = (i-R 2 ) 1/2 ^ + w, weL+, \\wf = R 2 <i. (201) 

To obtain estimate of E K ^ n from below note that (201) implies 

£o (T (*„)) = (1 - R 2 ) £o (*°) + So (w) > (1 - R 2 ) E° n + E° n \\w\\ 2 > E° n . 
Therefore, we obtain for L € {L}„, £ K>£ (L (B n )) < C 2n 

£ K , £ (T(iff ))>E°+ mm Q K ^ („) . 

*£E,£ reiS (u)<C 2 „ 

Inequality £ K ^(T(B n )) < C2n by Lemma 5.2 implies the boundedncss T (B n ) in 
H 1 . Using (iv) in Lemma 5.2 we conclude that 

E Kti , n > E° n - n 2 C 3n (1 + |ln K |) . (202) 

Combining (202), (198) and (199) we obtain (197). □ 

Now we prove the existence of solutions of eigenvalue problem with regularized 
logarithm. 

Theorem 5.7. Let — oo < £ < 0. Let integer n > 2 and k < 1 be so small that 

i - E°) + C ln n 2 + C 3n K 2 (1 + |ln K |) < 0, (203) 

where E® are defined by (191) and C\ n is the constant from (199). Then 
(i) the interval l n defined by (197) contains a critical point of 
(ii) there exists a solution iff n £ S ra d of (172) with eigenvalue u> = uj K ^_ n such 
that LE,^' S — lu k ^ n \ < ^- which satisfies 

Iwicf.n - Wo,n| < C n K 2 (1 + |ln/t|) . (204) 
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Proof. We follow here the lines of [9] with minor modifications. Point (ii) follows 
from point (i) and (186). To prove (i) assume the contrary, namely let n and k 
satisfy all the conditions and assume that that the interval /„ defined by (197) does 
not contain critical points of £ K ,£- We apply then Lemma 5.5 where we take 



1 

^2 



I - E° n ) + C lnK 2 + C 3nK 2 (1 + |ln«|) 

y 



with E' = E® — K 2 C3n (1 + |1hk|) and find a mapping rj : S ra d — > 3 ra d. By Lemma 
5.6 the value E K> £ >n defined by (196) or (200) lies inside /„■ By point (iv) of Lemma 
5.5 

r] ({u e H rad : £ K £ (u) < E K> £ >n + e}) c{ii£ S rad : (it) < E K ^ >n - e} . 

We can find T E {T} n such that max s „ P (*) < E K ^_ n + e/2. If we consider V (*) = 
r\ (r (^)) then we get maxg n T' (\I/) < E K ^ <n — e/2. To infer from this inequality a 
contradiction with definition (196) we have to check that V G {L}„ , namely to 
check (195). Note that on dB n 

£o((l-R 2 ) 1/2 K + 



and hence 



(1 - R 2 ) E° n + £ (v) <E° n + R 2 (E^ - E° n ) , 

£ K J(l-Rl) 1/2 *° 2 + v) (205) 



< E° + R 2 {El_ x E° n ) + g K ,t {(1 R 2 ) 1/2 K + 
£ K , e \ aB „ < E °n + R l ( E n-i - K) + (G K , e (*)) , (206) 



Therefore 



dB„ 

?2 



and using (199) and (203), since Rq = 1/9, we conclude that 

£ K ,t\ dBn < E° n + Rl {E° n _ x - E° n ) + C ln n 2 <E'-e. (207) 

implying r\ (^>) = ^ on dB n by point (i) of Lemma 5.5. Hence, V S {L}„ and we 
obtain a contradiction with (196). So we conclude that /„ contains a critical point 
Since £ K £ is diffcrentiable at vp', this critical point provides a solution to (172), 
and (204) follows from (180). □ 

Now we show that solutions of (172) decay exponentially as r = |x| — > oo. 

Lemma 5.8. Let n < no- Let < n < kq where kq is small enough for conditions 
of Theorem 5.7 to be fulfilled for n < hq. Let £ < satisfy 

l-y|e| + ^< 0. (208) 

Let P = § |£| ■ Let $f n € H ra d be a solution of (172) with eigenvalue ui — 
which is described in Theorem 5.7 with £ < £. Then there exist positive constants 
Ri and C which depend only on no, k and £, such that for £ < £, n < hq 

|*n Ml < C"e- pr /or r > Ri. (209) 

Proof. Since H^nll^-i are bounded uniformly in £ < and k < 1 we can use the 
following result of Strauss (see Lemma A. II in [9]) 

|*n Ml < C' a r- 1 ||tf n || H1 for r > R, (210) 
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where C and R are absolute constants, R > 1. We apply this inequality to radial 
solutions vf" = fy n of (172) with cj = uJ K ,£,n < and obtain 

|*«(r)| < Co^ 1 for r > R, (211) 

where, according to Lemma 5.2, the constant Co is bounded uniformly in k < 1, £ < 
and n. We write (172) in the form 

iv 2 * + V K ,^ = 0, (212) 



where 



^W = Ti-^ G '/^( |Vl/(x)l 



According to (211) using monotonicity of the function ln+,£ (s) with respect to s 
we obtain for |x| > R > R > 1, 

.2 



1 



(x) < ± + ■- (ln + ,, (^iT 2 C7 2 /C 2 ) + 3 



■ LJ. 



R ' 2 ^~ r * v " " ~ l " 9 
Note that 

(x) < I + y K + 3) + w if In ( K - 3 i?^ 2 C 2 /C 2 ) < £, R > R. 

Therefore, for any 0<«<Ko,£<0we can find R\> R such that 

ln( K - 3 i? r 2 C 2 /C 2 )<e, 
and if |£| is so large that (208) holds, we obtain for £ = £ 

k 2 

(x) < -— - \£\ +u < for |x| > R v (213) 

According to definition (160) V K> £ (x) is a monotonely increasing function of £, 
therefore this inequality holds for all £ < £. 

Note that \& (x) cannot change sign for |x| > R\. Indeed, if \& (x) = at 
|x| = R[ > Ri we could multiply (212) by \P and integrate over {|x| > R[} which 
would yield 

! fi|V*| 2 -^WI*| 2 )dx = 0, 

and (x) = for all |x| > R[ which is impossible for a radial solution of (172) 
with H'J'II = 1. So we assume (x) > for |x| > R\ (the case of negative ^ is 
similar). We compare the solution ^ of this equation in the domain |x| > R with 
the function 

Y = C Y -e- pr , p>0,C Y >0, 
r 

which is a solution of the equation 

V 2 F ~p 2 Y = 0. 
Subtracting half this equation from (212) we obtain 

iv 2 (tf - Y) + V kA (tf - Y) + \V K , e + ^ Y = 0. (214) 
We choose p 2 = ^ |£| and (213) implies 



K,« + \p 2 ) 
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We choose Cy > so that 

C Y e~ pRl > Co, 

where Co is the same as in (211). We assert that 

*(x) -Y (x) < for |x| > R 1 . (216) 

This is true at |x| = R\ according to the choice of Cy, the limit of ^ — Y as 
|x| — >• oo is zero. Assume now that (216) does not hold for some |x| > R\. Then 
^ (x) — Y (x) must have a point of local positive maximum at some |x| = R[ > Ri, 
at this point V 2 (* - Y) < 0, by (215) (V K>i + p 2 ) Y < and K,£ (* - Y) < 
by (213). This inequalities contradict (214), hence (216) holds. Inequality (209) 
follows from (216). □ 

Now we prove the main result of this section for the nonlinear eigenvalue problem 
with the original logarithmic inequality which corresponds to £ = — oo. 

Theorem 5.9. Let £ = — oo, uq > and let integer 2 < n < no and k < 1 be so 
small that (203) holds. Then there exists a solution \&„ £ S ra d n X of (172) with 
eigenvalue lu = uj K n with energy E* = uj K n + ^- which satisfies the inequality 

| < C n K 2 (1 + |ln k\) . (217) 

The solution ty n decays superexponentially as |x| — > oo, namely (209), where R\ 
depends on p, holds with arbitrary large p. 

Proof. We fix n and consider a sequence £ m — > — oo. We have a sequence of solu- 
tions ^ n ,m € S ra( j of (172) with , n satisfying (204) which exist according to 
Theorem 5.7 with £ = f m . Note that m according to (197) (#„, m ) € are 
bounded uniformly in m. Therefore ^ n ,m are bounded in H . Let us fix £ which 
satisfies (208) and take such m that £ m < £. All \t„ im satisfy (209). We can find 
a subsequence of ^ n ,m which converges to 'J/^oo G S ra d weakly in H , strongly 
in L 5 ({|x| < R}) n i 1 ({|x| < R}) for every R > and almost everywhere in K 3 ; 
we also choose the subsequence so that uj K ^ m ^ n — t uj Ki — oo,n G [£", £7"] - ^. We 
obtain that V 2 ^,,™ converges weakly to V vf n «, G S ra d. Now we show that ^„, !00 
is a solution of (172) with lj = u; Ki „ = w K — in the sense of distributions. The 
restriction of g K>im (* n>m (x)) to the ball {|x| < R} converges to g K -oo (* n ,oo (x)) 
in L 2 ({|x| < R}). Multiplying (172) with £ = £ m by an infinitely smooth test 
function / with compact support and integrating, we observe that we can pass to 
the limit in every term of the equation, obtaining that ^ n ,oo is a solution of (172). 
Now we show that \& raj00 G S ra( j n X. Since estimate (209) is uniform in £ m < £ and 
|G K ,-oo (C"e~ pr )| is integrable over R 3 we easily derive from dominated convergence 
theorem that 



Note also that 



as m — > oo. Indeed 



/ G re>Sm (*„ >m ) dx-> / (*„,„,) dx. (218) 

J|x|>_R! J|x|>i?! 

at 

/ G K , Cm (* n>m ) dx -> / G K ._oo (*„. oc ) dx (219) 

J|x|<i?i J\x\<Ri 



G K , U (|^„, m | 2 ) -GK,-oo(|*n 



G K , 4m (|*n, m | 2 ) - <?„,_«, (|*„, m | 2 ) + <?„,_«, (|*„, m | 2 ) - (\^n,oo\ 2 ] 
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Note that G K ,_oo \\^\ ) defined by (49) satisfies Lipschitz condition. 



Hence 



(l*n,J 2 ) 



<C||* n , TO |-|* n|00 || (l + | 



v& I 2 _ I* I 2 

^ n..m.\ * n. . rxn 



|x|<fl 



G' 



/t, — oo 1 I 1 rt.m 



G K 



^n,oo|| ( 1 ||^ f n,m||^4 



<C"||* 

^From (166) we obtain 

J\x\<R V 



n,oo 1 1 j^a 



(l^.ool 2 ) 

+ c" II* 



dx 



(220) 



n.m 1 n, oo r l 



G 



1 I ^n,m | 



dx < Ce«™ /4 



Li(|x|<7?.) 



I**,™. I dx. 



|x|<R 



(221) 

Using (220), (221) and convergence of *„, TO in L 1 ({|x| < R}) we obtain (219). 
From (218) and (219) we infer 



n,m 

)dx^ G K -oo (*„ j00 )dx. 
Jk 3 Jr 3 



(222) 



This implies that \&„ )00 G X with X defined by (137). From (222) and weak 
convergence in H 1 we infer that 

£ K (*n,oo) < lim SUp£ K (*„,m) • 
m— >-oo 

Since £ K ($„ im ) are bounded uniformly £ K (^n.oo) is bounded, hence ^ n ,oo G W. 
Therefore (172) holds in i?" 1 + X'. Similarly to (222) we obtain that ||*„, TO || 2 -> 
[1^71,00 1! 2 = 1, hence $ n>00 G 5 r ad- Since we can choose f arbitrary large, p in 
(209) can also be taken arbitrary large, and from convergence almost everywhere 
we derive that ^n.oo ( r ) decays superexponentially. □ 

The following theorem shows, roughly speaking, that there is no eigenvalues u> 
of the nonlinear problem in the gaps between small neighborhoods of eigenvalues of 
the linear problem.. 



Theorem 5.10. Let to satisfy the inequalities 



2 < w o,n + 5 <w < w , n+ i - S < 0, 



(223) 



where 5 > 0. There exists G4 (n) > sitc/i that if 

(5>C 4 (n)K 2 (l + |lnK|), (224) 

then the set of (radial or non-radial) solutions 3/ G 3 n X of equation (172) with 
such u is empty. 

Proof. Consider in L 2 (R 3 ) a linear operator (the classical Schrodinger operator 
with the Coulomb potential shifted by w ) 

1 9 1 

2 |x| 

Consider solution of (172) written in the form 



a,* = --<?«(*) 



(225) 
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with cj in the gap between two consecutive points of spectrum of O. Multiplying 
(172) (where £ = -oo) in H = L 2 (R 3 ) by # e 5 n X we obtain 

(O w * > *) = -i(ff« (*),*>, (226) 

where g K (iff) is defined in (181), (u, v) coincides with scalar product in H = L 2 (M 3 ) 
and determines duality between H 1 and H^ 1 and X and X' in (137). Note that 
9k {&) S -X 7 according to Lemma 9.3.3 in [15], hence the expressions are well-defined. 
Since 

( ^- |*| 2 dx < C e ||*|| 2 + e||*|| H i||*|| for * G H 1 
J |x| 

with arbitrary small e > 0, we have 

(o«*,*> > -c£||#|| 2 - 

Using (128) we estimate (g K (iff) , iff) and we obtain for solutions of (172) the follow- 
ing inequality: 



In / |V*| 2 rfx 



\\\n 2 m-c"\\n 2 <l^ 

Hence||*|| ff i < M for solutions >t6SnX 

Let E® C H 1 (R 3 ) be the finite-dimensional space with the orthonormal basis 
of eigenfunctions $fj of the operator Oq = O with negative eigenvalues u)j < wo >n 
(since we are not restricted to radial functions, the eigenvalues u>j < Wj+i may have 
high multiplicity and cb n < wo,n for n > 2). Let If n be orthoprojection in L 2 (M 3 ) 
onto i?^. We have by (223) 

((o u -o u n n )%-<si) > Kn+i-^)ll*ll 2 . 

Since the eigenfunctions of O arc bounded, continuous and decay exponentially, 
we see that operators 

n n 

n„* = J2 (% CM = ]T (d&j - W ) (*, ^ 

are bounded from L p (R 3 ) to L p (R 3 ) for any given p,p', oo > p > 1, oo > p' > 1. 
Solutions $ £ S (1 1' are bounded in if 1 , and by Sobolev imbedding in i 6 (R 3 ) . 

Since | 5 „ (<f>)\ < Ck 2 (1 + |1hk|) (\<f>\ 3/2 + |*| 1/2 ) we obtain 

||n nflK (*)|| <C" K 2 (l + |ln K |) (227) 
Multiplying (172) by * - E„* we obtain 

{O u V, (* - En*)) = -i (<?« (*) , (* - !!„*)) 



Therefore 



(w 0l „+i - w) || (1 - E„) *|| 2 < -I ( 5k (*) , (tf - n n *)) (228) 



<~\ j 9, W*dx+i||E r ,g K (*)||||*| 



Using (128) and (227) to estimate the right-hand side we obtain 

(wo,n+i - w) || (1 - E n ) *|| 2 < C lK 2 (1 + |ln K |) . (229) 
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Multiplying (225) by U n ^ we obtain 

K„- W |||n„*|| 2 < |(M,n n *)| = i|(n B ^(*),n n *>| (230) 

< ^\\n n g K (m <C 2K 2 (l + |ln K |). 
^From (229) and (230) using (223) we obtain 

£||*|| 2 < C 3 « 2 (1 + |M) 
If we take in (224) C4 (n) = 2C3 we obtain a contradiction with ||\I>|| = 1. Hence, a 
solution of (172) with such ui does not exist. □ 

Remark 6. Theorem 5.9 states that for every energy level of the linear Schrddingcr 
hydrogen operator there exists an energy level of nonlinear hydrogen equation with 
the Coulomb potential with an estimate of difference of order ~ k 2 In k. Based on 
(173) we see that if x = h where h is the Planck constant and n = then uj 1 = 
coincide with standard expression for hydrogen frequencies, and a natural choice 
for parameter x is X = ft- 

Remark 7. A theorem similar to Theorem 5.9 (without supcrexponcntial decay 
of solutions) can be proven not only for the logarithmic nonlinearity, but for more 
general subcritical nonlinearities which involve a small parameter k. Moreover, 
if corresponding functions g K are of class C 1 and their derivatives satisfy natu- 
ral growth conditions, then the standard bifurcation analysis shows that every E® 
generates a branch E K>n of solutions & n £ S ra( j of (172) for small K. 

Remark 8. In Theorem 2.1 we assumed in the macroscopic case electron size a to 
be very small and in this section we assume k = a±/a to be very small, which seem 
to contradict one another. But a closer look shows that in fact the both cases are 
consistent one with another, if we assume the possibility a\ — > 0. To show that 
this assumption is quite reasonable, we look at the values of physical constants. 
Note that a typical macroscopic scale is a macr ~ 10 _3 m and the Bohr radius a\ ~ 
5.3 x 10 -11 m. Note also that the error in (82) is of order a 2 /a 2 nacr <C 1 and here 
k 2 = a\/a 2 <C 1 which is possible if a\ <C a 2 <C a^ cr Since af/a^acr ~ 2.8 x 10~ 15 
we consider a 2 /a 2 nacr — > as a reasonable approximation. 

Remark 9. In the case when the energy functional with the Coulomb potential 
£cb ( V I / ) is considered not only on radial functions, the eigenvalues of the linear 
problem are not simple and have very high degeneracy. Note that if we start from the 
relativistic equation (13) and look for w-harmonic solutions we arrive to a nonlinear 
perturbation of so-called relativistic Schrddinger equation, see [34] . The spectrum of 
this linear problem has a fine structure and the degeneracy with respect to the total 
angular momentum disappears. If we fix the value of z-component L z of the angular 
momentum L, for example by looking for solutions of the form ^1 =e 1T ™ a ^xo ( r , z ) 
with a fixed integer m z , one can eliminate the remaining degeneracy in the linear 
problem. This provides a ground to expect that the lower energy levels of the 
nonlinear relativistic Schrddingcr equation can be approximated by the lower energy 
levels of the linear relativistic Schrddinger equation if n is small. 

Remark 10. Note that the decoupled equation (152) with b = for proton has a 
unique ground state solution (gausson) which corresponds to the minimum E\ of 
energy (see Theorem 3.4 ) and has a sequence of solutions with higher critical energy 
levels E?, j > 2 (see [14]). Since the ground state is unique, E? > Ef if j > 1. 
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The logarithmic nonlinearity G' 2 = — K P m y K P 3 l"0| 2 /Cg J — 3 for proton 

involves the parameter k p which is different from the parameter k for electron. A 
simple rescaling and use of (122) shows that dependence of E% (n p ) on k p is very 
simple: E? (k p ) = K p E^ (1). There are physical motivations based on experimental 
estimates of the size of the proton to assume that n p is large, k p 3> 1. So we can 
assume that the gap E p - (n p ) — Ef (k p ) > {E? (1) — Ef (1)) n p is larger than 1 for 
j > 2. (Numerically calculated values of E% (1) — E\ (1) are given in [13], they are 
of order 1). At the same time ujQ. n given by (191) satisfy inequality — ^ < <Jo,« < 0. 
Consequently, since the total energy £ defined by (144) of the system (141)-(143) 
equals the sum of energies of electron and proton one may expect for small b and k 
and large n p that the lower negative energy levels of the functional £ cannot branch 
as b varies from E? + ujq.u with j > 2 and must branch from E P + uJq,„. Hence we 
expect that if 6 and K are small and n p is large, the lower negative energy levels 
of the functional £ which corresponds to the system (141)-(143) are close to lower 
discrete energy levels of £cb which are described in Theorem 5.9. 
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